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SUMMARY 


Stress intensity levels at various parts of the NASA B-52 carrier aircraft pylon were examined for the case 
when the pylon store was the space shuttle solid rocket booster drop test vehicle. Eight critical stress points were 
selected for the pylon fatigue analysis. Using fracture mechanics and the half-cycle theory (directly or indirectly) 
for the calculations of fatigue-crack growth, the remaining fatigue life (number of flights left) was estimated for each 
critical part. It was found that the two rear hooks had relatively short fatigue life and that the front hook had the 
shortest fatigue life of all the parts analyzed. The rest of the pylon parts were found to be noncritical because of their 
extremely long fatigue life associated with the low operational stress levels. 


INTRODUCTION 

The NASA B-52 carrier aircraft has been used to carry various types of test vehicles for high-altitude drop tests. 
Test vehicles carried by the B-52 aircraft in the past include X-15, HL-10, highly maneuverable aircraft technology 
(HiMAT), and drones for aerodynamic and structural testing aircraft (DAST), and the space shuttle solid rocket 
booster drop test vehicle (SRB-DTV). The test vehicle was attached to the B-52 pylon (airborne launch system) 
by one front hook and two rear hooks (see fig. 1, in which the store is the SRB-DTV). The B-52 pylon (fig. 2) 
is almost 30 years old and has been subjected to considerable structural fatigue. Recently, the two rear hooks of 
the pylon fractured during the towing of the B-52 aircraft carrying the SRB-DTV on a relatively smooth taxiway. 
During towing, the hook loads were far below the design limit loads. Careful observation of the fracture surfaces 
of the two failed rear hooks revealed that “micro” surface cracks existed at the rounded inner boundaries of the two 
failed rear hooks (Ko and Schuster, 1985). This incident raised the serious concern about the remaining fatigue life 
of the front hook and critical parts of the pylon and also the available fatigue life of the new rear hooks (made of 
high-performance material and of better design). Before resuming the new series of SRB-DTV drop tests, it was 
necessary to conduct proof tests of the pylon by loading it up to specified limit loads to establish confidence in 
the pylon structural integrity. The main purpose of proof tests is as follows. During the proof tests, if some pylon 
components should fail, then those components are to be replaced. If all the pylon parts survive the proof tests, then 
fracture mechanics can be applied to estimate the fatigue life (or number of flights left) for each of the critical parts 
of the pylon based on proof test data and data obtained from the SRB-DTV test flights following the proof tests. 
This report shows the application of fracture mechanics and the half-cycle method to calculate the remaining fatigue 
life (number of flights left) for each critical component of the pylon carrying the SRB-DTV as a store. 

NOMENCLATURE 


CCLAMP 

CC LH 

EOS 

EPSL 

FLIDAB 

HMS 

NASTRAN 

PCM 

SMAX 


compute crack length at maximum points 

computed crack lengths 

end of segment 

epsilon 

flight database 

subroutine; h, min, sec 

NASA structural analysis 

pulse code modulation 

sigma maximum 



SMIN sigma minimum 

SRB-DTV space shuttle solid rocket booster drop test vehicle 

SYMBOLS 


A 

Ab 


C 

c 

Dcl 

Dcr 

Dj 

E 

Ft 

f 

9 

Ip 

Ir 

Iy 

i 

i 

Ki 

Fl c 

Kj 

Umax 

k 

i 

Mk 

My 

M 2 

m 

N 


crack location parameter 

pylon-wing attachment bolt cross-sectional area 

crack length of edge crack, half crack length of through-thickness crack, or depth 
of surface crack 

initial fictitious crack length (or depth) established by the proof tests 

initial fictitious crack depth established by the proof tests for the pylon-wing 
attachment bolts 

limit crack length (or depth) associated with the operational load level 

crack size (length or depth) after the £th flight of random stress cyclings 

material constant in Walker crack growth rate equation 

half-length of surface crack 

left drag load 

right drag load 

pylon-wing attachment point drag loads 
complete elliptic function of the second kind 
number of flights left after the Mi flight 
fraction of limit stress (or limit load) 
gravitational acceleration 
SRB-DTV pitching moment of inertia 
SRB-DTV rolling moment of inertia 
SRB-DTV yawing moment of inertia, Iy ft* Ip 

integer associated with half-cycles, or associated with different regions of crack growth curve 

integer associated with half cycles 

Mode I stress intensity factor 

Mode I critical stress intensity factor 

function of hook loads and pylon accelerations 

Mode I stress intensity factor associated with cr max , K max , = AM K^maxs/^o/Q 
modulus of elliptic function 
integer associated with flights 
flaw magnification factor 

aerodynamic pitching moment acting on SRB-DTV, in.-lb, positive nose up 
aerodynamic yawing moment acting on SRB-DTV, in.-lb, positive nose right 
Walker exponent associated with K max 
number of stress cycles available for operations 
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Py 

Pz 

V 
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9 
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P<li » Rbi , Rc 1 
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Sa 

ScL 

Sj 
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Sn, > S pi , Sq { , 
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T 

t 

V A 

Vbl 

Vbr 

Vcl 
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V; 

W 

x,y,z 

X 

o; 

OiRUS 

0 

A 

Sdi 


A K 

Vx 


number of stress cycles used during £th flight 
Walker exponent associated with stress ratio R 
aerodynamic drag force acting on SRB-DTV, lb, positive aft 
aerodynamic sideforce acting on SRB-DTV, lb, positive right 
aerodynamic lift force acting on SRB-DTV, lb, positive up 
roll inertia load factor, p = —6 X / 9 
surface flaw shape and plasticity factor 
pitch inertia load factor, q = —6 y /g 
dynamic pressure 

coefficients of hook load equation for hook load i? t - 

stress ratio, R — Gmin/ ® max 
yaw inertia load factor, f = —9 z /g 
front side load 
rear side load 

pylon-wing attachment point side loads 
coefficients of stress equation for cr, 


time 

plate thickness 

front hook vertical load 

left rear hook vertical load 

right rear hook vertical load 

left drag shaft vertical load 

right drag shaft vertical load 

pylon-wing attachment point vertical loads 

weight of SRB-DTV 

rectangular Cartesian coordinates 

location of aerodynamic center of pressure on x axis for SRB-DTV 
slope of crack growth curve (or crack growth rate) 
fuselage angle of attack 
yaw angle 

amount of crack growth during the l\h flight 

crack growth increment caused by one cycle of constant-amplitude stress 
cycling under loading magnitude of A Ki and R x 

Mode I stress intensity range, A K = AMk( cr TOQX — a mtn ) , s/ixa/Q 
= —x/g, inertia load factor in the i direction at SRB-DTV center of gravity 



Vy 

Vz 

Vx 

Vy 

Vz 

Ox 

Oy 

O z 
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O-oo 

<^Vj 

CTy 

O’D' 

^mai 


CTi 


[ ]t 
[•],["] 


= —y/g, inertia load factor in the y direction at SRB-DTV center of gravity 

= —z/g, inertia load factor in the 2 direction at SRB-DTV center of gravity 

dimensionless acceleration in the 1 direction at pylon center of gravity 

dimensionless acceleration in the y direction at pylon center of gravity 

dimensionless acceleration in the 2 direction at pylon center of gravity 

angular displacement about the x axis (roll angle) 

angular displacement about the y axis (pitch angle) 

angular displacement about the 2 axis (yaw angle) 

angular coordinate of critical points <74 and cr 6 

remote uniaxial tensile stress 

uniaxial tensile stress induced by the proof load 

uniaxial tensile stress associated with the operational load level 

uniaxial tensile stress at pylon-wing attachment point induced by the proof loads 

yield stress 

tensile strength 

maximum stress of the stress cycle 
minimum stress of the stress cycle 
tensile stress at stress point i 
limit value of cr, established by proof tests 
ultimate shear strength 

angular coordinate for semielliptical surface crack 
quantity associated with the £th SRB-DTV test flight 
time derivatives 


FRACTURE MECHANICS 


The upper part of figure 3 shows the most common types of cracks: through-thickness crack, surface crack, and 
edge crack. According to fracture mechanics, the stress intensity factor Ki for the Mode I deformation (tension 
mode) associated with any type of crack can be expressed as 


Ki = AMk<Joo 



( 1 ) 


where 


A 

is 

M k 

is 

CToo 

is 

a 

is 


the crack location parameter (A = 1 for the through-thickness crack, A = 1.12 for 
both the surface and the edge cracks (see fig. 3)), 

the flaw magnification factor (for a very shallow crack, Mr = 1 ; as the depth of the 
crack reaches the back surface of the plate, Mr = 1 .6 (see fig. 3)), 

remote uniaxial tensile stress, 

half crack length for the through-thickness crack, or length of the edge crack, or 
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depth of the surface crack (see fig. 3), and 
Q is surface flaw shape and plasticity factor, 


i 2 


Q = [E(k)] 2 -0.212 — 

Lay 

where cry is the yield stress and E( k ) is the complete elliptic function of the second kind defined as 

rir/2 j 

E(k)=J y 1 — fc 2 sin 2 <f> d<f) 

where the angle (f> is defined in figure 3 and the modulus k of the elliptic function is defined by 

2 TO] 2 

k 2 = 1 - 1 


( 2 ) 

(3) 


(4) 


where c is the half-length of the surface crack (see top center of fig. 3). The lower part of figure 3 shows the plots of 
a/2 c as a function of Q for different values of ffoo/ cry and the plot of Mk as a function of a/t, where t is the plate 
thickness. 


PROOF LOAD TESTS 

The purpose of proof load tests is to load the entire aircraft structure (or its components) up to certain proof load 
levels to test structural integrity and to establish initial fictitious crack sizes for critical structural components for 
fatigue life analysis. The proof load levels are usually slightly lower than the design limit load conditions associated 
with different maneuvers. If a previously undetected crack in a certain component exists and is larger than the critical 
crack size, that component will certainly fail during the proof load tests and will be replaced. Thus, a catastrophic 
accident during flight could be avoided. If the entire structure survives the proof load tests, then the critical point of 
the structural component has been subjected to a proof stress induced by the proof loads. If Ki c is the critical 
stress intensity factor (or material fracture toughness) of the structural component material, the maximum crack 
length af the structural component can carry under the proof loads without failure (or rapid crack extension) may 
be calculated from equation (1) by setting Ki = Ki e and <Joo = cr£>. In reality, cracks may not develop during the 
proof tests; however, it will be assumed that a fictitious crack of length af has been created at the critical point 
of the structural component during the proof load tests. During actual operations, the structural component will be 
subjected to much lower stress levels than the proof stress level cr£, . If a^ is the peak operational stress level (highest 
peak of the stress cycles), then according to equation (1) the structural component can carry a fictitious crack of size 
a°, which is much larger than af . The value a° thus determined will be considered as the limit crack size toward 
which the initial crack af is allowed to grow after repeated operations. The crack size difference, a° - af , will then 
be the available crack size increase permitted for the structural component in repeated operations. The left side of 
figure 4 shows the plot of crack length, a, as a function of normalized stress CToo/°ft, where ay is the failure stress 
of the material. The lower the operational stress level, the more the available crack size increase is allowed for the 
structural component. The right side of figure 4 shows the corresponding plot of the crack length, a, as a function 
of number of constant stress cycles N. 

REMAINING FLIGHTS 

If the structural component is cycled under constant stress amplitude (idealized case for the purpose of discussion) 
for N\ cycles during the first flight with the associated crack growth of Aai , then the number of remaining flights 
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F\ may be estimated approximately from (see fig. 4) 


F\ 



Aoi N\ 


( 5 ) 


Equation (5) is based on the assumption that the crack growths for all subsequent flights will be equal to A oi of the 
first flight. In reality, the crack growth rate will steadily increase with the number of flights accumulated, and the 
actual remaining flights wifi be less than the value F\ predicted by equation (5) if the number of remaining flights 
is large (that is, F\ 1 ). For a moderate range of Fi , equation (5) should give a sufficiently accurate prediction 
of the number of remaining flights. The amount of crack growth Aoi in equation (5) may be calculated from the 
following Walker equation (Forman and others (1967), Walker (1970), and Schijve (1979)) for fatigue-crack growth 
rate under constant-amplitude stress cyclings: 


-777 = C(fr moi ) m (l -f?) n = C(A20 m (l -i?)"-"* (6) 

aN 

where C, m, and n are material constants and K max , A K, and R are respectively the stress intensity factors associ- 
ated with c max t the stress intensity range, and the stress ratio, given by 


Kmax = AM K Omax^m = A K( 1 - R) ~ l (7) 

na 

AK = AKd fci&rnax 0>mn) (8) 

(9) 

'O max 

where o ma x and o m in are respectively the maximum and minimum stresses of the constant-amplitude stress cycles. 

In reality, the stress cycles at the critical points of the structural components are not constant-amplitude stress 
cycles. To apply equations (6) to (9) to the variable-amplitude stress cycles, certain methods must be developed. In 
this report, a half-cycle theory (Starkey and Marco (1957), Incarbone (1963), and Smith (1965)) is used to calculate 
the fatigue-crack growth rate for variable-amplitude stress cycles. 


HALF-CYCLE THEORY 

The upper part of figure 5 shows an example of random stress cycles (variable-amplitude loading history). The 
stress history curve is the combination of a series of increasing and decreasing load half-cycles of different loading 
magnitude (A K, R), as shown in the lower part of figure 5. The half-cycle (or half-wave) theory (Starkey and 
Marco (1957), Incarbone (1963), and Smith (1965)) states that the damage (or crack growth) caused by each half- 
cycle of either increasing or decreasing load is assumed to be equal to half the damage caused by a complete cycle of 
the same loading magnitude (A K, R). This means that the damage caused by the complete cycle could be equally 
divided between the two phases of increasing and decreasing loads. The loading sequence can thus be resolved into 
half-cycle groups of increasing and decreasing loads (lower part of fig. 5). Each half-cycle (either increasing or 
decreasing load) can then be considered as a half-cycle of the constant-amplitude cyclings under the same loading 
magnitude (A K, R) and can be computed separately in time sequence to estimate the corresponding damage. The 
half-cycle theory thus permits accurate evaluation of the load spectmm from a recorded load time history. If a* (Z = 
1, 2, 3, . . . ) is the final crack length after Ng,(Z= 1, 2, 3, ... ) random stress cyclings in the £th flight, then according 
to the half-cycle theory, at may be calculated from 

a= a*_, + Jj — = a £ _! + Aa* (£=1,2,3,...) (10) 

i=i 1 
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with 


2 Ni c , 

Aae = ^2 ~T~ ^ a<_i = of (11) 

i=i 2 

where 6 a, /2 is the crack growth increment induced by the ith half-cycle under the loading magnitude A if,- and i?,, 
and which is assumed to be equal to the crack growth increment induced by a half-cycle of the constant-amplitude 
stress cycling fatigue test under the same loading magnitude ( A Ki, Ri ) . Thus, by using equations ( 6 ) to (9), 6o,/2 
may be calculated from 


6o< 

~ 


2 


" da ' 
_dN_ 


= J [( Kmax ) »] ”* ( 1 - ft) n = |(Air,) m (l - H,) n-Tn 


( 12 ) 


The crack length, o, associated with A Ki (see eq. (7)) will be the summation of initial crack size of and all the 
crack growth increments created by all the previous half-cycles. Namely, 


a 


%- 1 


;=i z 


(i> 2 ) 


(13) 


Similar to the constant-amplitude stress cyclings (see eq. (5)), the number of remaining flights Fi predicted using 
the random stress cyclings of the first flight will be 



(14) 


Similarly, after the £th flight, the number of remaining flights Fg calculated from the flight data of the £th flight 
will be 


Fg=^ 


- Q£-l 
A a i 




1 


IS Sat 

h 2 


ag - at - 1 
a i — a^_i 


(15) 


Figure 6 graphically illustrates how to evaluate the crack growth increment 6o,/2 associated with the ith half-cycle 
of the random stress cyclings by using the plots of A Ki as a function of da/dN for different values of R. 


STRESS POINTS 


From the stress analyses of the B-52 pylon hooks (Ko and Schuster, 1985) and pylon parts (North American 
Aviation, Inc., 1963), eight stress points (stresses cti to os) were selected for the pylon fatigue analysis. Each of the 
stress points is located at a critical point of each of the eight parts. Table 1 lists the names of the eight pylon parts and 
the associated material and fatigue properties of the parts. The stress points < 71 , 02 , and <13 are respectively located 
at the critical points of the inner circular boundaries of the front hook, left rear hook, and right rear hook (figs. 7 
and 8 ). The stress point <74 is located at the critical point of the boundary of the left drag shaft (fig. 9). The stress 
point <75 is located at the critical point of the right drag shaft (fig. 10 ). The stress points 0* and 07 are respectively 
located at the critical points of the left drag pin (fig. 9) and right drag pin (fig. 10). Finally, the stress point as is 
located at the lower edge (fig. 11 ) of the right drag shaft lower support. 

Based on the NASA structural analysis computer program (NASTRAN) stress analyses of the three B-52 pylon 
hooks (ai , <72 , and <73 ; (Ko and Schuster, 1985)), the proof tests of the left drag shaft (<74 ) and the stress analyses of 
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the B-52 pylon parts (<75 to <7g; (North American Aviation Inc., 1963)), the eight stresses <71 to ag (in kips, or 10 3 
lb/in. 2 ) can be expressed in terms of the hook loads (in pounds) as follows: 


<ri = 7.3522 x 10~ 3 Vk 


02 = 5.8442 x 10 -3 Vsl 

03 = 5.8442 x 10 — 3 Vbr 

04 = (6.3481 -5.4514 x 10 _5 Dcl) x 10 ~ 3 D C l cos 9 

-(9.9653+ 3.4363 x 10 ~ 4 S C i) x lO" 3 S C Lsin0 


where the angular coordinate 9 of critical point is given by 


9 = tan 1 


ScL 


Dei 


= tan 


-1 


8210 


= 14 .36 ‘ 


proof 32080 

05 = 2.1868 x 10" 3 D C r 

06 = 7.1420 x 10 ~ 3 (Dcl cos 9 - Scl sin 9) 

07 = 7.3055 x 10 ~ 3 D C r 
og =4.2141 x 10 ~ 3 D C r 


(16) 

(17) 

(18) 

(19) 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where the front hook vertical load Va, the left rear hook vertical load Vbi,. the right rear hook vertical load Vbr, 
the rear side load Scl > the left drag load Dcl, and the right drag load Dcr, may be (1) obtained from strain gage 
measured data if the strain gage instrumentations are available for measuring the hook loads or (2) calculated from the 
hook load equations given in appendix A using accelerometers and aerodynamic data if the strain gage data are not 
available. The stress points at the pylon-wing attachment bolts and several other pylon parts were also investigated; 
however, those stress points turned out to be noncritical, and therefore they are not included in the pylon fatigue 
analysis. Appendix B shows the calculations of the pylon-wing attachment bolt tensile loads and the initial surface 
crack sizes associated with the attachment bolt stress points. 


CALCULATIONS OF CRACK SIZES 


Table 2 shows different combinations of the maximum hook loads used in the proof tests. To obtain the highest 
possible values of the stresses a, (t = 1, 2 , .... 8) for the fatigue analysis, the highest values of the maximum hook 
loads or the highest values of the combined applied loads (for 04 and <76 only) were used. The values of <7 ,• (i = 1 , 
2 , .... 8) thus obtained were used to calculate the initial fictitious crack size af and the operational limit crack size 
a° for each stress point using the following equations (eq. (1) rewritten): 



a 


O _ 

c 


Q 

7 T 


Ku 


AM k o* 

Kic 

AMfcfa* 


(»= 1 , 2 ,.. ., 8 ) 


7 2 


(i = 1 , 2 ,..., 8) 


(24) 

(25) 


where / is the fraction of the limit stress a* (that is, fa* is the peak operational stress) and the critical stress intensity 
factors Ki c for all the stress points are given in table 1. 

In the calculations of af and a° (see eqs. (24) and (25)), all the surface cracks will be considered to be semi- 
elliptical in shape, with aspect ratio a/2 c = 1 /4 . This value is based on the worst case of the failed old rear hooks. 
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For the value a/2 c = 1/4 , the lower left plot of figure 3 gives the surface flaw shape factor of Q = 1 .25 taking 
<?<x>/vy - 1 because most of the limit stresses a* (i = 1 , 2 ,..., 8) calculated from equations (16) to (23) exceeded 
the corresponding yield points of the parts materials. Only for the stress point ag , Q was taken to be Q = 1 because 
the crack was a two-dimensional edge crack. Also, the value of the flaw magnification factor Mr was set to unity 
( Mk - 1) for all the stress points because the depths of the cracks were very small compared with the characteristic 
dimensions (such as thickness, diameter, and depth of the pylon parts under consideration). Table 3 shows the 
calculated initial fictitious crack size of established from the proof tests, the amount of crack growth per flight, and 
the number of flights left calculated from the £th flight data for the eight stress points. As mentioned previously, the 
limit crack size a° for the actual operation is based on the peak stress recorded during the flight. Since the operational 
peak stress varies each flight, the number of flights left calculated from the £th flight data may be larger or smaller 
than those predicted for the earlier flights (such as cri and crj). 

FLIGHT DATA 

A total of eight successful (that is, with SRB-DTV launching) SRB-DTV flight tests were conducted. One 
additional flight, which occurred between flights 4 and 6, was an aborted flight without the launching of the SRB- 
DTV. During the entire SRB-DTV flight-test program, there were only three channels available for recording the 
strain gage outputs. For flights 1 to 7 (including the aborted flight), the three available channels were used to measure 
V A , Vbl> and Vbr ( Dcl> Dcr . and Scl were not measured). At the end of flight 7, the left drag shaft was removed 
from the B-52 pylon for proof test and for the installation of strain gages for measuring Dcl and Scl • For flights 
8 and 9, the two channels used for measuring Vbl and Vbr were switched to measure Dcl and Scl ( Vbl and 
Vbr were not measured). The value of Dcr was not measured during the flight-test program. Accelerometer data 
(accelerations in x, y, and z directions) were available for the entire test program. However, the accelerometer 
outputs were poor and unreliable and therefore could not be used to calculate accurate hook loads. An alternative 
method (called the indirect method, using strain gage measured hook loads from other flights) for calculating the 
crack growths without using accelerometer data is described later in this report. 

CALCULATIONS OF LOAD SPECTRA 

To perform the crack growth analysis, the loading histories for the eight stress points must first be calculated. 
Using the strain gage measured values of Va, Vbl> and Vbr (for flight 1) and Dcl and Scl (for flight 9), the eight 
stresses cr, (i = 1,2,..., 8) may be calculated from equations (16) to (23). Figures 12 to 19 show portions of the 
loading histories calculated for the eight stress points during taxiing for flight 1 (cri to cr?) and for flight 9 (a a to 
as ). These stress history curves (load spectra) were filtered down to 5 Hz to eliminate the small-amplitude, high- 
frequency cycles, which are not important in the present fatigue analysis. Notice from figures 12 to 19 that the load 
spectra for all the stress points are random. It must be mentioned that in the calculations of <75 , at, and <78 (see 
figs. 16, 18, and 19 and eqs. (20), (22), and (23)), the strain gage measured values of Dcl for flight 9 were used 
instead of Dcr values which were not measured during the entire SRB-DTV flight-test program. Since the values of 
Dcl and Dcr are quite close and their perturbation parts are mutually out of phase by a half cycle (see eqs. (A-23) 
and (A-24) in app. A), using the Dcl strain gage data for calculating crs , < 77 , and erg should give reasonably accurate 
results for the crack growth analysis. 

Appendix C discusses in detail the computer programs used to calculate <r t - (i = 1, 2, . . . , 8) from the strain 
gage measured hook loads or from the accelerometer output (if the data are good) using the hook loads equations 
presented in appendix A. 
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CRACK GROWTH CALCULATIONS 


Direct Method 

For certain flights, if the load spectra for certain stress points can be calculated direcdy from the strain gage 
measured hook loads, the half-cycle method may be used “directly” to calculate the crack growths for these stress 
points. By resolving the calculated load spectra (see figs. 12 to 14) into a series of half-cycles with different loading 
magnitudes (A K,R) and applying the half-cycle theory mentioned previously, the crack growth increments per half- 
cycle could be calculated in time sequence and summed to give the total crack growth per flight for each stress point. 
Finally, equation (14) or (15) may be used to calculate the number of remaining flights associated with each stress 
point. It must be emphasized that in using the half-cycle theory, every half-cycle in the load spectrum is calculated 
without neglecting insignificant low-amplitude half-cycles. Thus, the half-cycle method gives accurate evaluations 
of the crack growth. Based on the status of the flight data (strain gage measured hook loads) mentioned previously, 
the direct method for calculating crack growth may be used for the eight stress points in the following way: 


< 7 i flight 1 to flight 9 

< 72 , <73 flight 1 to flight 7 

<74 to <78 flights 8 and 9 (after the left drag shaft proof tests and the strain gage instrumentation) 


Indirect Method 

The indirect method for calculating the crack growths for several stress points must be used for some flights: 


<72, <73 flights 8 and 9 

<74 to <78 flight 1 to flight 7 


The crack growth curve for each stress point for a typical flight may be divided into four major regions with 
different slopes (crack growth rates) and different time durations: 

1. taxiing, 

2. takeoff run, 

3. airborne, cruise power, and 

4. cruising. 

The upper curve of figure 20 shows the crack growth curve calculated for, say, <74 using the strain gage data of 
flight 9 and using of as the initial crack size (the initial crack size for flight 1). From this curve, the four averaged 
slopes <*,• (i = 1 , 2, 3, 4) for the four regions can be determined. Using the <*, thus obtained and the known durations 
Ti (» = 1, 2, 3, 4) for the four regions of flight 1, it is possible to construct a crack growth curve for flight 1, Gower 
curve of fig. 20). Thus, the estimated crack growth at stress point <7,- for flight 1 may be calculated from 

4 

A ai = ][>,?;■ (26) 

t=i 

For calculating the amount of crack growth at stress point <74 for flight 2 , one has to go back to the strain gage data 
of flight 9 , recalculate the crack growth curve, and determine the new set of slopes a,. This time, the initial crack 
size of of + A 01 must be used. This process can be used up to flight 7 . Calculations of crack growths for the stress 
points <75 to as may be carried out in a similar manner. In applying the indirect method for the calculations of crack 
growths at stress points <72 and <73 for flights 8 and 9 , the master slopes or,- (crack growth rates) established from the 
strain gage data of flight 7 were used. 
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Missing Data Extrapolations 

During the flights, short intervals of strain gage data were lost as a result of noise, interference, and the in- 
terruption by parts of the B-52 aircraft in the line-of-sight signal path from the transmitter on the payload to the 
ground receiver. To account for those missing data in a particular maneuver region (such as taxiing, takeoff run, 
airborne and cruise power, or cruising), “forward” extrapolation of the crack growth curve with respect to time was 
performed using the average crack growth rate a, (or slope) from the previous maneuver transition point (or check 
point). Figure 21 illustrates this missing data extrapolation. 

For certain flights, the initial taxiing data were not recovered. For such a case, the “backward” extrapolation 
of the crack growth curve had to be performed. As shown in the upper part of figure 22, the crack growth curve 
was first calculated starting from point B to point B' (the next check point, or maneuver transition point) using the 
proof crack size af. Then the averaged slope ai of the crack growth curve lying between points B and B' was 
established. Using the slope a\ thus established, the amount of crack growth A a that occurred during time T\ of the 
initial missing data region could be calculated from 


Ao = aiT) (27) 

After A a was calculated, point B was then shifted upward to location B" by an amount A a (that is, the crack 
size for point B is increased to af + A a; see lower part of fig. 22), and then the final calculations of the crack growth 
curve proceeded from point B" . 

Appendix C discusses in detail the computer programs used in the calculation of fatigue-crack growths using the 
half-cycle method directly and indirectly. 

RESULTS 

Figures 23 to 30 show the crack growth curves calculated for the eight stress points cr* (» = 1,2,..., 8) for flight 1 . 
The crack growth curves for other flights are quite similar and therefore are not shown. Figures 23 to 25 (curves for 
stress points cri to CT3 ) were calculated by using the direct method; figures 26 to 30 were constructed by using the 
indirect method. The maneuver transition points are indicated in those figures. In figures 23 to 25, the data dropout 
regions are also shown. For the front and the two rear hooks (at , <J2 , <73 , figs. 23 to 25) the crack growth rate is 
most rapid during the initial stage of taxiing and the takeoff run and becomes very slow during cruising because of 
relatively low-amplitude stress cyclings. However, for the rest of the stress points (04 to ag , figs. 26 to 30) the crack 
growth rate is very low during taxiing because of negligible drag and becomes larger during takeoff and cruising 
because of increased drag. Table 3 summarizes the results of the crack growth analysis of the eight pylon parts 
calculated for each SRB-DTV test flight. For practical applications, the calculated values of Ft should be reduced 
by using a safety factor in the range of 2 to 4 to account for uncertainties. Notice in table 3 that the value of Ft does 
not decrease monotonically with accumulated flights. The reason for this is that the duration and flight conditions 
for each flight are different (different A at), and therefore the value of Ft predicted at the end of the £th flight could 
be smaller or larger than the value Ft - 1 predicted at the end of (t - 1 )th flight, depending on the size of crack growth 
A at during the t\h flight. Figure 31 shows plots of the number of flights left F\ as functions of percent limit loads 
for the eight stress points for the first flight. The arrows in the figure indicate the locations of the peak operational 
load levels. Notice that the left and right drag shafts (0-4 , as), the left and right drag pins (ae, a-j ), and the right 
drag shaft support (ag) have long fatigue life left after the completion of the SRB-DTV flight tests because of low 
operational stress levels. These parts (a4 to ag) turned out to be noncritical because the operational crack sizes a° 
exceeded the characteristic dimensions of the corresponding parts (see table 3). Fatigue lives for the three hooks 
(at, a2, a3) after the final flight are relatively short. The front hook (ai) has the shortest fatigue life because of 
relatively high operational stress levels, and it is therefore the most critical part of the pylon. Figure 32 shows plots 
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of amounts of crack growth with the increase in (light number for the front hook (the most critical part) and the two 
rear hooks (critical parts) during the entire SRB-DTV test program. 

CONCLUSION 

Eight stress points, located respectively at eight critical points of the B -52 pylon parts, were selected for the pylon 
fatigue life analysis. Using fracture mechanics and the half-cycle theory (directly and indirectly) for the calculation 
of fatigue-crack growth, the remaining fatigue life (or number of flights left) was estimated for each of the pylon 
critical parts. The results showed that the left and right drag shafts (stress points <74 . 0 ’s), and drag pins (stress points 
<76 , <77), and the right drag shaft lower support (stress point as), have long fatigue life because of low operational 
stress levels and can be considered as noncritical parts. Also, the fatigue lives of the three hooks (stress points <71 , 
( 72 , <73 ) are relatively short. The front hook (stress point crj ) has the shortest fatigue life compared with other pylon 
parts because of relatively high operational stress level, and it is therefore the most critical part of the pylon. 


National Aeronautics and Space Administration 

Ames Research Center 

Dryden Flight Research Facility 

Edwards, California 

July 25, 1986 
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Table 1. Material properties for critical parts (stress points). 


Stress 

point 

Part 

name 

Part 

number 

Material 


j 

ksi 

cry, 

ksi 

0-1 

Front hook 

2581-630620 

Inconel 718 alloy 

175 

145 

O ' 2 

New left rear hook 


— 

AMAX MP35N alloy 

250 

235 

CJ 3 

New right rear hook 


— 

AMAX MP35N alloy 

250 

235 

0*4 

Left drag shaft 

2581-630706 

4340 steel 


180 to 200 

225 

0-5 

Right drag shaft 

2581-630707 

4340 steel 


160 to 180 

163 

0*6 

Left drag pin 

2581-630718 

AMAX MP35N alloy 

250 

235 

07 

Right drag pin 

2581-630716 

AMAX MP35N alloy 

250 

235 

0*8 

Right drag shaft 

2581-630711 

4340 steel 


160 to 180 

163 


lower support 








Stress 

Part 

r, 


c, 



Ki c , 


point 

name 

ksi 

^(ksiV ^)-" 1 

m 

71 

ksi \/im 

<y\ 

Front hook 

135 

9.220 

x 10 -1Z 

3.60 

2.16 

125 


0*2 

New left rear hook 

141 

2.944 

x 10~ u 

3.24 

1.69 

124 


0*3 

New right rear hook 

141 

2.944 

x 10- 11 

3.24 

1.69 

124 


0*4 

Left drag shaft 

156 

1.043 

x 10 " 9 

2.70 

2.27 

105 


0 'S 

Right drag shaft 

108 

1.043 

x 10 “ 9 

2.70 

2.27 

105 


0-6 

Left drag pin 

141 

2.944 

x 10" u 

3.24 

1.69 

124 


0-7 

Right drag pin 

141 

2.944 

x 10" n 

3.24 

1.69 

124 


0*8 

Right hand drag 

108 

1.043 

x 10“ 9 

2.70 

2.27 

105 



lower support 









Table 2. Proof test hook loads. 

Maneuver 

identification 

number 




Hook load, lb 



V A 

Vbl 

Vbr 

Sa 

Sd 

Dcl 

Dcr 

1 

35,019 

28,708 

31,437 

1,684 

5,369 

25,660 

26,113 

2 

35,603 

36,754 

36,088 

0 

0 

0 

0 

3 

27,416 

32,859 

42,756 

3,995 

-10,099 

0 

0 

4 

27,415 

42,871 

32,744 

3,995 

10,099 

0 

0 

5 

29,728 

31,806 

30,991 

0 

-8,210 

32,080 

30,060 




Maximum proof load 




36,520 

44,110 

44,230 

3,590 

10,100 

31,900 

30,030 


-9,650 -9,660 
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Table 3. Summary of crack growth analysis. 


Stress 

point 


Proof 

load 

test 

49.48 

Flight duration, min 
64.70 67.55 52.83 

99.43 

1 

2 

Flight number 
3 

4 

5 

CTl 

cti 

0.0990 

0.1008295 

0.1021277 

0.1044369 

0.1062016 

0.1120766 


< 


0.3333052 

0.3333052 

0.2835354 

0.2835354 

0.2059647 


A at 


0.0018295 

0.0012982 

0.0023092 

0.0017646 

0.0058250 


<T°c-<Tl 


0.2324758 

0.2309243 

0.1790985 

0.1773338 

0.0938881 


F t 


128 

179 

78 

101 

17 


Percent limit load 


54.50 

51.31 

59.09 

55.77 

69.33 

&2 

CTl 

0.0734 

0.0739887 

0.0746227 

0.0753477 

0.0760176 

0.0776691 


< 


0.2076090 

0.2057364 

0.2057364 

0.1856398 

0.1776967 


Aa* 


0.0005887 

0.0006340 

0.0007250 

0.0006699 

0.0016515 


a.° - at 


0.1336203 

0.1311137 

0.1302389 

0.1096222 

0.1000276 


Ft 


227 

207 

180 

164 

61 


Percent limit load 


59.46 

59.73 

59.47 

62.88 

64.27 


CTl 

0.0730 

0.0737705 

0.0742656 

0.0751271 

0.0758267 

0.0775514 


a° c 


0.2037274 

0.2037274 

0.1930696 

0.1930696 

0.1930696 


A at 


0.0007705 

0.0004951 

0.0008615 

0.0006996 

0.0017247 


a° - at 


0.1299569 

0.1294618 

0.1179425 

0.1174803 

0.1155182 


Ft 


169 

263 

137 

168 

78 


Percent limit load 


59.86 

54.19 

61.49 

60.68 

58.79 

<J 4 

CTl 

0.1410 

0.1412347 

0.1415575 

0.1416992 

0.1419603 

0.1424520 


at 


3.4865426 

3.4865426 

3.4865426 

3.4865426 

3.4865426 


A a* 


0.0002347 

0.0003228 

0.0001417 

0.0002610 

0.0004917 


a° c - a* 


3.3453079 

3.3449851 

3.3448434 

3.3445823 

3.3440906 


Fl 


14,255 

10,362 

23,608 

12,813 

6,802 


Percent limit load 


20.11 

20.11 

20.11 

20.11 

20.11 

0’S 

CTl 

0.3815 

0.3816163 

0.3817765 

0.3818368 

0.3819657 

0.3822081 


at 


14.0298566 

14.0298566 

14.0298566 

14.0298566 

14.0298566 


A ai 


0.0001163 

0.0001602 

0.0000604 

0.0001289 

0.0002424 


(T° c -<T t 


13.6482403 

13.6480801 

13.6480198 

13.6478909 

13.6476485 


Ft 


117,418 

85,211 

226,097 

105,904 

56,317 


Percent limit load 


16.49 

16.49 

16.49 

16.49 

16.49 

0-6 

CTl 

0.0872 

0.0872515 

0.0873233 

0.0873520 

0.0874094 

0.0875180 


at 


3.7446171 

3.7446171 

3.7446171 

3.7446171 

3.7446171 


A a* 


0.0000515 

0.0000719 

0.0000286 

0.0000574 

0.0001086 


a° c -ai 


3.6573656 

3.6572938 

3.6572651 

3.6572077 

3.6570991 


Ft 


70,992 

50,872 

127,739 

63,628 

33,669 


Percent limit load 


15.26 

15.26 

15.26 

15.26 

15.26 

CTJ 

at 

0.1010 

0.1010611 

0.1011467 

0.1011772 

0.1012458 

0.1013746 


at 


3.7143264 

3.7143264 

3.7143264 

3.7143264 

3.7143264 


A at 


0.0000611 

0.0000855 

0.0000305 

0.0000683 

0.0001291 


a° c ~ at 


3.6132653 

3.6131797 

3.6131492 

3.6130806 

3.6129518 


Ft 


59,105 

42,259 

118,442 

52,942 

27,980 


Percent limit load 


16.49 

16.49 

16.49 

16.49 

16.49 

erg 

at 

0.1750 

0.1751166 

0.1752773 

0.1753379 

0.1754673 

0.1757108 


at 


6.4201312 

6.4201312 

6.4201312 

6.4201312 

6.4201312 


A at 


0.0001166 

0.0001607 

0.0000606 

0.0001294 

0.0002435 


a° c -at 


6.2450146 

6.2448539 

6.2447933 

6.2446639 

6.2444204 


Ft 


53,609 

38,884 

103,134 

48,282 

25,660 


Percent limit load 


16.51 

16.51 

16.51 

16.51 

16.51 



Table 3. Concluded. 


Stress 

point 


Proof 

load 

test 

95.55 

Flight duration, min 
58.55 44.23 

67.65 

6 

Flight number 
7 8 

9 

O-l 

a* 

0.0990 

0.1154196 

0.1186624 

0.1207648 

0.1267761 




0.2059647 

0.2059647 

0.2059647 

0.2059647 


A at 


0.0033930 

0.0032429 

0.0021024 

0.0060133 


a° c - at 


0.0905455 

0.0873023 

0.0851999 

0.0791886 


Ft 


27 

27 

41 

14 


Percent limit load 


54.09 

55.42 

53.49 

66.72 


at 

0.0734 

0.0790510 

0.0799232 

0.0806427 

0.0816830 


a °c 


0.1776967 

0.1440195 

0.1440195 

0.1440195 


A Ql 


0.0013819 

0.0008722 

0.0007195 

0.0010403 


a? - a t 


0.0986457 

0.0640983 

0.0633768 

0.0623365 


Ft 


72 

74 

89 

60 


Percent limit load 


63.12 

71.39 

71.39 

71.39 

<73 

at 

0.0730 

0.0790016 

0.0799302 

0.0807175 

0.0818549 


a °c 


0.1930696 

0.1779445 

0.1779445 

0.1779445 


A at 


0.0014502 

0.0009286 

0.0007873 

0.0011374 


a° c - at 


0.1140680 

0.0980143 

0.0972270 

0.0960896 


Ft 


79 

106 

124 

85 


Percent limit load 


59.73 

64.05 

64.05 

64.05 

(74 

a t 

0.1410 

0.1428617 

0.1431225 

0.1412063 

0.1415535 


a °c 


3.4865426 

3.4865426 

3.4865426 

3.4865426 


A at 


0.0004097 

0.0002609 

0.0002063 

0.0003473 


a c — at 


3.3436809 

3.3434201 

3.3453363 

3.3449891 


Ft 


8,162 

12,817 

16,220 

9,633 


Percent limit load 


20.11 

20.11 

20.11 

20.11 

<75 

at 

0.3815 

0.3824089 

0.3825361 

0.3826368 

0.3828103 


a? 


14.0298566 

14.0298566 

14.0298566 

14.0298566 


A at 


0.0002008 

0.0001272 

0.0001007 

0.0001735 


a°-at 


13.6474477 

13.6473205 

13.647220 

13.6470463 


Ft 


67,990 

107,295 

135,614 

78,668 


Percent limit load 


16.49 

16.49 

16.49 

16.49 

O '. 6 

at 

0.0872 

0.0876037 

0.0876603 

0.0872441 

0.0873218 


< 


3.7446171 

3.7446171 

3.7446171 

3.7446171 


A at 


0.0000857 

0.0000566 

0.0000441 

0.0000778 


a°- a t 


3.6570134 

3.6569568 

3.6569127 

3.6568349 


Ft 


42,634 

64,571 

82,938 

47,002 


Percent limit load 


15.26 

15.26 

15.26 

15.26 

<77 

at 

0.1010 

0.1014763 

0.1015435 

0.1015959 

0.1016893 


°c 


3.7143264 

3.7143264 

3.7143264 

3.7143264 


A at 


0.0001018 

0.0000672 

0.0000524 

0.0000934 


a°- at 


3.6128501 

3.6127829 

3.6127305 

3.6126371 


Ft 


35,508 

53,779 

69,005 

38,684 


Percent limit load 


16.49 

16.49 

16.49 

16.49 

<7g 

at 

0.1750 

0.1759127 

0.1760407 

0.1761420 

0.1763168 




6.4201312 

6.4201312 

6.4201312 

6.4201312 


A at 


0.0002019 

0.0001280 

0.0001013 

0.0001748 


a° - at 


6.2442185 

6.2440905 

6.2439892 

6.2438144 


Ft 


30,948 

48,797 

61,675 

35,742 


Percent limit load 


16.51 

16.51 

16.51 

16.51 



Forward 



Figure 1. Geometry of space shuttle solid rocket booster drop test vehicle (SRB-DTV) attached to B-52 pylon. View, looking 
inboard at right side of B-52 and SRB-DTV. 
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Figure 2. Geometry of B-52 pylon. 
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Variable-amplitude 
stress cycles 



Figure 5. Resolution of random stress cycles into half stress cycles of different stress ranges. 
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Figure 7. Front hook and the location of critical stress 
point cri . 



Figure 8 . Rear hook (left or right) and the location of critical 
stress point or a 3 ). 
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Figure 9. Left hand drag shaft and drag pin and the location of critical stress points 0-4 and <X6 . 
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Figure 12. Stress cycles for critical stress point <ti at front hook V A . 
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Figure 13. Stress cycles for critical stress point CT 2 at left rear hook Va l 
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Figure 20. Construction of crack growth curve for flight 1 using slopes calculated from strain-gage data of flight 9. 
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Figure 23. Crack growth curve for stress point 1 a \ , flight 1 (direct method). 
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Figure 24. Crack growth curve for stress point 2 cj 2 , flight 1 (direct method). 
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Figure 29. Crack growth curve for stress point 7 a-j , flight 1 (indirect method). 



Figure 30. Crack growth curve for stress point 8 erg , flight 1 (indirect method). 


Crack 

length, 

in 


Crack 

length, 

in 


35 





APPENDIX A 

EQUATIONS FOR B-52 PYLON HOOK LOADS 

The following equations for the B-52 pylon hook loads (forces acting on hooks) were derived in accordance with 
the dimensions and the reference system shown in figures A- 1 to A-3. All the inertia forces and moments are written 
in terms of inertia load factors at the SRB-DTV center of gravity (rj x , r\ y , rj z , p, q, r ), respectively, where 


rh = 




P = 




(A-l) 

(A-2) 


1. Front hook vertical load Va (positive downward) 

Considering the moment about line CL - CR (see fig. A-l), 

Va = Jifo {9 ‘ 85( ~ WTlx + Px) + 60 8 W( 1 + Vz) + lfp9+W(60.8 2 + 9 .85 2 )]q - P 2 x - M y } ( A-3) 

where W is the SRB-DTV weight; P x aerodynamic drag force, and P z aerodynamic lift force acting on the 
SRB-DTV; Ip SRB-DTV pitching moment of inertia; g gravitational acceleration; x location of aerodynamic 
center of pressure on x axis for SRB-DTV; and M y aerodynamic pitching moment on SRB-DTV. 

2. Front side load Sa (positive inboard) 

Considering the moment about point 0 (see fig. A-2), 

Sa = 22^29 {6 °' 8 Wr ly-^y9+W(60.8) 2 ]r-P y x-M z } (A-4) 

where 220.9 is the effective moment arm (defined in North American Aviation, Inc., 1963); Iy the SRB-DTV 
yawing moment of inertia; P y , aerodynamic sideforce acting on the SRB-DTV; and M z aerodynamic yawing 
moment acting on the SRB-DTV. 


3. Left rear hook vertical load Vbi (positive downward) 

Considering the moment about the line passing through point BR and parallel to the x axis (see fig. A-3), 

Vbl= 52 ( 26 - 43 1 + P*) -Va-Pz] -9. 85(Wp y - P y ) -16.3125^ 

- [Ir9 + W( 26 .437 2 + 12.163 2 )]p) (A-5) 


where Ip is the SRB-DTV rolling moment of inertia. 

4. Right rear hook vertical load Vpp (positive downward) 
Considering the forces in the z direction (see fig. A-3), 


Vbr = 1 + Vz ) - Vbl — V A — Pz 


( A-6) 


5. Rear side load Scl (positive inboard) 

Considering the forces in the y direction (see fig. A-2), 


Scl = Wr\ y - P y - Sa 


(A-7) 



6. Left drag load Dcl (positive rearward) 

Considering the moment about the line passing through point 0 and parallel to the z axis (see fig. A-2), 

D CL = mT475 {60 - SWriy ~ [lY9+ M7(60 - 8)2]f - p ^- M -> + j(-Wvx+Px) (A-8) 
where the effective moment arm 1484.475 is defined in North American Aviation, Inc. (1963). 

7. Right drag load Dqr (positive rearward) 

Considering the forces in the x direction (see fig. A-2), 

D C r = - Wt) x + P x - Dcl ( A-9) 

The aerodynamic forces and moments appearing in equations (A-3) to (A-9) may be calculated from the 
following equations (Quade (1977), configuration 1): 

P x = 5.019840 x 10 3 g (A- 10) 

P y = (1.59264 -0.521165/3) x 10 3 g (A-ll) 

P z = ( —0 .220320 + 0 .378778 cxfus) x 10 3 g (A- 12) 

M y = -(3.385053 + 0.477254 cx FUS ) x 10 5 q (A-13) 

M z = (-2.085160 + 0.636736/3) x 10 5 g (A-14) 

where q is dynamic pressure (lb/in. 2 ), cxfus fuselage angle of attack, and /3 yaw angle. 

For the SRB-DTV, 

W = 48,143 lb (A- 15) 

I P g = 1.4591 x 10 9 lb-in. 2 % I Y g (A-16) 

I R g = 4 .5537 x 10 7 lb-in. 2 (A- 17) 

Using these numerical values, the hook load equations may be written as 


V A = — 2247 .434 t\ x + 13872.485(1 + tj a ) + 7.780750 x lO 6 ? 

+ 0 .046682 P x - 0 .288152 P z - 0 .004739 M y (A- 18) 

S A = 13287 .261 7] y - 7.431309 x 10 6 f - 0.275996 P y 

-0. 004539 M z (A- 19) 

V B l = 1 123. 717 r} x - 13067 .992 r? y + 17135 ,258( 1 + tj,) 

- 1.632632 x 10 6 p- 3.890375 x 10 6 <j + 2 .292681 x 10 6 f 

- 0 .023341 P x + 0 .271448 P y - 0 .355924 P z 

+ 0 .002370 My — 0 .001400 M z (A-20) 

V BR = 1123 .717 Tj x + 13067 .992 rj y + 17135 ,258( 1 + tj 2 ) 

+ 1.632632 x 10 6 p- 3.890375 x 10 6 g - 2 .292681 x 10 6 f 

- 0 .023341 P x - 0 .271448 P y - 0 .355924 P 2 

+ 0.002370 M y + 0.001400 M z (A-21) 

S C l = 34855 .739 rj y - 7.431309 x 10 6 f - 0.724004 P y 

+0 .004539 M z (A-22) 

D cl = -24071.50Tb + 1971. 804 7j y - 1.102792 x 10 6 f 

+ 0 .5 P x - 0 .040957 P y - 0 .000674 M z (A-23) 

-24071 .50 ri, - 1971. 804 rj y + 1.102792 x 10 6 f 

+ 0 .5 P x + 0 .040957 P y + 0 .000674 M z (A-24) 
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igure A-l. Location of SRB-DTV and B-52 pylon front and rear hooks. Dimensions in inches. 
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APPENDIX B 

PYLON LOADS AT B-52 WING ATTACHMENT POINTS 
AND THE ASSOCIATED STRESS POINTS 

The pylon is attached through bolts to the B-52 wing at four points (designated 1 to 4) as shown in figures 2 and 
B-l. The pylon vertical loads Vj (j = 1 to 4), pylon side loads Sj, (j ' - 1 to 4), and pylon drag loads D 3 and £? 4 
are expressed in terms of the hook loads and the pylon accelerations. After the pylon loads are calculated in terms 
of the proof load combinations given in table 2, the stress levels at the four attachment points are examined to see 
whether those stress points are critical in the pylon fatigue analysis. 

For Case 1 of Quade (1977) (that is, 1/3 of AKs is carried by bolts 1 and 2), the pylon loads may be written as 


Vi = -0 .53115150 K\ - 0 .05045848 X 2 - 0.18451742 X3 

+ 0.00118936X4 -0.01006607X5 - 0.00363813 X 6 (B-l) 

V 2 = -0 .53115150X1 -0.0564432 X2 + 0.21055295 X3 - 0.00152446 X 4 

+ 0 .01462529 X 5 - 0 .00363813 X 6 (B-2) 

V 3 = 0 .02946747 Xi + 0 .01234619 X 2 - 0 .02494134 X 3 

+ 0 .00163321 X 4 - 0 .02899156 X 5 + 0 .00344146 X 6 (B-3) 

V 4 = 0 .03233554 Xi + 0 .05375659 X 2 - 0 .00109420 X 3 

- 0 .00129812 X 4 + 0 .02443235 X 5 + 0 .00383481 X 6 (B-4) 

S\ = —0 .0468 1422 X 2 - 0 .53 124658 X 3 + 0 .00364924 X 4 (B-5) 

s 2 = -0 .04681422 X 2 - 0 .53124658 X 3 + 0 .00364924 X 4 (B-6) 

S3 = 0 .04681422 X 2 + 0 .03124658 X 3 - 0 .00364924 X 4 (B-7) 

Sa = 0 .04681422 X 2 + 0 .03124658 X3 - 0 .00364924 X 4 

Dj — Da — -0 .500000 X 2 (B-8) 

where X ; 0 = 1 to 6) are functions of hook loads and pylon accelerations: 

Ki=V a + Vbl + Vbr + V CL + V C r - 11 70 77, (B-9) 

X 2 = D C l + D C r - 1170% (B-10) 

X 3 ==Sa + S C l+ 1170% (B- 11) 

X 4 = -78 .5 S A + 132 .5 Sex + 1170% x 53.239 - 1170% x 13.5 

+ 44 . 812 D C r- 17. 812 D C l (B-12) 

X 5 = 25.32475 S A + 1170% x 9.13725 + 41. 63725 S C l 

+ 26 ,437( Vs£, - V B r) + 31.312 (V C l - V CR ) (B-l 3) 

X 6 = 78 .5 V A + H70% x 53 .239 - 132 ,5( V BL + V B r + Vcl + Vcr. ) 

+ 1170% x 13.793 -46.293 (D C l + D C r) (B-14) 


where Vcl and Vcr are respectively the left and right drag shaft vertical loads and %, %, and % are respectively 
the nondimensional pylon center-of-gravity accelerations in the x, y, and z directions (see fig. B-l). 

Substituting equations (B-9) to (B-14) into equations (B-l) to (B-4), the four pylon vertical loads Vj O' = 1 to 4) 
may be expressed as 

V! = —0 .816745 ^4 - 0.5328035^ — 0. 315216 Vbl + 0. 096777 D C l 

- 0 .446051 Scl + 0 .217017 V B r + 0 .171259 D C r 

- 18.460970% - 249.412739%+ 394.829484% (B-l 5) 
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Vi = — 0 .816745 Va + 0.700604875^ + 0.33754952 7^ + 0.13913043 D C l 
+ 0 .61751886 S C l - 0 .43574806 V B r + 31 .4059356 % - 4 .96397070 % 

+ 394 .829487 % - 0 .46432619 D C r (B- 1 6) 

7 3 = 0 .29962208 7,4 - 0 .88735234 5,4 - 1 .19297585 V B i + 0 .33992389 7 BB 

- 1. 015669845 Scl -0.17605610 Dcl - 0.07378191 £> C h + 15.29603010% 

- 237 .3859800% + 179 .8935723 % (B-17) 

7 4 = 0.3338681374 + 0.719551385^ + 0. 17064125 7 fl £- 1 .12119483 7 bb 

+ 0. 84420077 5cl - 0.10064616 - 0.1 8 193931 Do* 

- 21.5135856% + 179.0564256% + 200.4513147% (B-18) 

Now, based on various combinations of the proof hook loads given in table 2, the magnitudes of Vj (j = 1 to 4) 
induced by the proof loads will be examined. For the proof tests, % = % = % = 0 since the proof tests are static 
in nature. 


Maneuver 1. 

— For maneuver 1, table 2 gives 

74 = 35,019 lb , V BL = 29,708 lb , V BR = 31,437 lb 

5,4 = 16841b, S C l = 5369 lb 
D C l = 25,660 lb , D C r = 26,113 lb 
Using those values, equations (B-15) to (B-18) yield 


7i = -27,165 lb 

(B-19) 

Vi = -24,544 lb 

(B-20) 

Vj = -21,943 lb 

(B-21) 

7» = -20,246 lb 

(B-22) 


Maneuver 2. 

— For maneuver 2, table 2 gives 

74 = 35,603 lb , V BL = 36,754 lb , V BR = 36,088 lb 


Sa = Scl = Dcl = Dcr = 0 
Using those values, equations (B-15) to (B-18) yield 

7i = -32,832 lb (maximum for 7i ) (B-23) 

Vi = -32,398 lb (B-24) 

7 3 = -20,912 lb (B-25) 

7» = -22,303 lb (B-26) 
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Maneuver 3. 


— For maneuver 3, table 2 gives 

V A = 27,416 lb , V BL = 32,859 lb , V BR = 42,756 lb 

S A = 3995 lb , Scl = -10,099 lb , D C l = D C r = 0 
Using those values, equations (B-15) to (B-18) yield 

Vx = -21,0951b 
V 2 = -33,3691b 
V 3 = -9739 lb 

V4 = —38,828 lb (maximum for V4) 


(B-27) 

(B-28) 

(B-29) 

(B-30) 


Maneuver 4. 

— For maneuver 4, table 2 gives 

V A = 27,415 lb , V BL = 42,871 lb , V BR = 32,744 lb 


S A = 3995 lb , Scl = 10,099 lb 
Dcl = Den = 0 

Using those values, equations (B-15) to (B-18) yield 

V\ - —30,973 lb (B-31) 

Vi = -13,7121b (B-32) 

V 3 = —45,602 lb (maximum for V 3 ) (B-33) 

V 4 = -8844 lb (B-34) 


Maneuver 5. 

— For maneuver 5, table 2 gives 

V A = 29,728 lb , V BL = 31,806 lb , V BR = 30,991 lb 

S A = 0 , Scl — -82101b 
Dcl = 32,080 lb , D CR = 30,060 lb 
Using those values, equations (B-15) to (B-18) yield 


Vx = -20,229 lb 

(B-35) 

V 2 = —41 ,634 lb (maximum for Vi ) 

(B-36) 

U 3 = -18,0291b 

(B-37) 

U 4 = -35,023 lb 

(B-38) 


In the preceding equations, the negative signs in front of V) 0 = 1 to 4) denote that the pylon reaction forces V) 
are pointing downward, giving tension to the attachment bolts. 
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Stress Points 


The four stress points oyi , ay 2 . &V 3 , and cry 4 are selected at the surfaces of four attachment bolts 1, 2, 3, and 
4, respectively (fig. B-l). In the calculations of the initial fictitious surface crack sizes at the stress points, the peak 
absolute values of Vj (j = 1 to 4) will be used. They are (see eqs. (B-23), (B-36), (B-33), and (B-30)) 


Vi = 32,832 lb (B-39) 

Vi = 41,634 lb (B-40) 

V 3 = 45,602 lb (B-41) 

V 4 = 38,828 lb (B-42) 


The four attachment bolts each have a diameter of 1 in., and the cross-sectional area of each attachment bolt is 
Ab = tt( 0 .5) 2 = 0 .785 in. 2 Thus, the stresses ovi to <7v4 may be calculated as 


ovi = y- = = 41 .824 ksi (B-43) 

Aq 0.785 

ai ' 2= ^ = ^ =53o37ksi ^ 

<t ” = ^ = ^ =58 0921:51 <b - 45) 

(B - 46) 

The four attachment bolts are NAS636AH 90 type bolts made of 4340 steel and heat treated for 1 80 to 200-ksi tensile 
strength. From table 1 , the critical stress intensity factor Ki for the four attachment bolts will be K / = 105 ksi \ZirT 
Taking Q = 1 .45 (since 0 < CToo/ cr y < 0 .6 ), Mk = l.A = 1.12, the initial fictitious surface crack sizes , a p v2 . 
a^ 3 , and ay 4 associated with the four stress points oyi , ayi . OV 3 , and ay* may be calculated from equation (1) as 


p _ Q ( if/ \ 2 _ 1-45 / 105 

Wl t IT Vl. 12 X 1 x 41.824 

p _Q ( Ki \ 2 _ 1.45 / 105 

V2 ‘k\AMk<% 2 ) * Vl- 12 x 1 x 53.037 

p _ Q f Ki \ 2 _ 1.45 / 105 

173 7 t\AM k o p vz ) 7 r V 1. 12 x 1 x 58.092 

p _Q ( if/ \ 2 _ 1.45 / 105 

' 74 TT\AM K <y p VA ) 7 r Vl- 12 x 1 x 49.462 


2 .3 19 in. 
1 .442 in. 
1 .202 in. 
1 .658 in. 


(B-47) 

(B-48) 

(B-49) 

(B-50) 


Those crack sizes are much greater than the diameter of the attachment bolts (1 in.) because the stress levels <jy X 
to are relatively low. Thus, the four attachment bolts are not critical and will not be used in the pylon fa- 
tigue analysis. 
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v' i ew looking foevard 

Figure B-l. Rear view of B-52 pylon showing hook loads and loads at pylon-wing attachment points. Dimensions in inches. 



APPENDIX C 

COMPUTER PROGRAMS USED IN THE CALCULATIONS OF a, 

AND CRACK GROWTH 


Crack Growth Calculation Overview 

Background. 

A method of testing various vehicles at NASA is to carry the test specimen aloft under the wing of a B-52 
carrier aircraft and release it to descend while transmitting measurements to ground based receivers. One such item 
is the SRB-DTV, which is used to develop improved parachute designs. While under the B-52 wing the specimen 
is supported by three hooks (one forward and two aft) and several sway braces and drag pins, all parts of a pylon 
structure. To determine the life expectancy of these parts, the stresses at the critical locations of each are computed 
from the accelerometers and strain gages on the payload. These stresses are then used to compute a theoretical crack 
length increase during taxi, takeoff, and flight up to the point of payload release. 


Problems. 

Certain problems must be taken into account when calculating crack growth. 

1. Physical placement of the strain gages and the accelerometers is on the test specimen, not on the pylon, and 
they are separated by some distance and through flexible structures from the stress sites on the pylon. Rotation 
and translation calculations are therefore required to compute the stress at the desired location. 

2. Vibration, especially in the accelerometer data, creates many false peaks that must be removed by heavy 
filtering. 

3. Data loss occurs from various causes, but the greatest losses are believed to be caused by interference by parts 
of the B-52 aircraft when they come between the antenna on the payload and the ground based radio receivers. 
One loss occurs at the beginning of most flights between engine start and when the aircraft has taxied to a 
point out of the “radio shadow” of the rooftop antenna. Losses of 0.1 sec or more are compensated for by 
extrapolation. 

4. Data file size is such that a full flight cannot be processed in the computer at one time. Flight durations are ap- 
proximately 45 min (or more), the sampling rate is 1000 samples/sec and 68 parameters are recorded, creating 
more than 183 million computer words, not counting ancillary information such as parameter descriptions, 
calibrations, and tables of dropouts. Processing the flight tape requires breaking the flight into several files (a 
minimum of two and a maximum of seven have been used to date). 


Procedure. 

General purpose programs: A series of general-purpose programs use calibration, configuration, and setup 

files to convert the raw pulse code modulation (PCM) data into a flight database (FLIDAB) file in engineering units 
(fig. C- 1). All these standard processes are shown as one in figure C-2 in the block labeled PROGDAB PROGRAMS. 
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Special purpose programs: Five special-purpose programs were specifically written to calculate and display 

crack growth. Because of the developmental nature of the programming requirements, some of the special programs 
are loaded from cards and compiled for each pass. 

1 . SRBHOOK uses the general-puipose DABIR subroutine to extract the required data channels and then applies 
translation, rotation, and moment arm transformations and low-pass filtering to determine the stress at each of 
the desired locations. Input to SRBHOOK controls the part of the flight file attached and is used to eliminate 
large areas of data losses. Each of the attached segments is put into a separate file, identified as SIGMA files 
in figure C-2. 

2. The THIN program attaches the SIGMA files one at a time and recatalogs them after reducing the number of 
data points and interpolating missing data sections. 

3. CONCAT attaches all the THIN files for a particular flight and concatenates them into a single file. 

4. CCLAMP then performs the crack length calculations and computes the number of flights remaining on each 
part. Two files are created: One file contains the crack length for each site at specified time intervals during 
the flight, the other file is a list of start and stop times for each segment of the flight where data were filled in 
by extrapolation. These are identified as CRACK LENGTH and DROPOUTLIST files in figure C-2. 

5. Plots are created by CPLOT, which utilizes the CRACK LENGTH and DROPOUTLIST files to show the 
progression of the cracks at periodic intervals during the flight and to identify extrapolated sections. 

! 

I 

File naming convention: A seven-character file name is used in which each character provides significant 

information. A typical file name, S0705T3, would indicate S series, flight 07 (first since modification of critical 
parts), 0 suffix created by the THIN program, and third segment of the flight. 

1 . The first character indicates a computation series. All data for this report were generated in the S series. 

2. The second and third characters constitute a two-digit flight-number code. The flight numbers used in this 

I report are obtained by subtracting 6 from this number (the total number since this series of flights started). 

3. The fourth character is zero for most flights but is used as a suffix when two operations have the same flight 
number (for example, flights 11 and 11 A). 

4. The fifth character indicates the filter cutoff frequency to identify files produced as a result of experimentation 
with other frequencies. 

5. The sixth and seventh characters identify the program that generated the file, as in the following tabulation: 


Sn SRBHOOK 
Tn THIN 
CX CONCAT 
CL CCLAMP cracklengths 
DL CCLAMP dropout list 


6. The seventh character is used by SRBHOOK and THIN to identify the section of the flight represented (that 
is, 1 denote the first part, 2 the second part, . . . , A the 10th, B the 1 1th, and so on). 

Program running time: The following estimated program running times on the CDC Cyber 170-730 (Control 
Data Corp., Minneapolis, Minnesota), for each minute of flight are given in the following: 
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Program 

Min CPU time 

Min of flight 

SRBHOOK 

7.75 

1 

THIN 

1.5 

1 

CONCAT 

0.1 

1 

CCLAMP 

3.5 

1 

CPLOT 

a 



“Running time is more dependent on the amount 
of annotation than on flight time. CPU time has 
varied between 30 and 50 sec. 


SRBHOOK Program 

SRBHOOK is a program used to read nine parameters from an SRB flight data base (FLIDAB), to filter these 
data, and then to compute several stresses <x,-, which are stored on a permanent file. This file is used to compute hook 
crack growth through the THIN and CCLAMP programs. The nine parameters read from the FLIDAB are 

1 . MOl , X -axis acceleration in g; 

2. M02, Y -axis acceleration in g\ 

3. M03, Z-axis acceleration in g; 

4. M04, roll rate in deg/sec; 

5. M05, pitch rate in deg/sec; 

6. M06, yaw rate in deg/sec; 

7. M56, forward hook load in kips; 

8. M67, aft hook load in kips; and 

9. M69, aft hook load in kips. 

Card input to SRBHOOK is in the form of three name lists: 

1 . NAMELIST /SRBNAM/ NSIGMA, BREAK, 1ST, IET 

NSIGMA is the number of stresses ct, to be computed (eight were computed so far), 

BREAK is break frequency (in Hertz) used with the fourth-order low-pass Butterworth filter used in 
SRBHOOK (a value of 5 was usually used), 

1ST is requested FLIDAB start time (hr, min, sec, msec), and 

IET is requested FLIDAB end time (hr,min,sec,msec). 

2. NAMELIST /TABLE1/ RA,RB,RC,RD,RE,RF 

These are seven-element arrays of coefficients used to compute seven hook loads i?, as shown in table C-l. 
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3. NAMELIST /TABLE2/ SK,SL,SM,SN,SP,SQ,SR 

These are eight-element arrays of coefficients used to compute eight stresses <x,- as shown in table C-2. 

The SRBHOOK program begins by reading the three name lists. Flight data procession system subroutines 
DABOP and DABIR are then called to extract the nine FLIDAB parameter values at each time point of the 
desired interval. The data are extracted from the FLIDAB at a rate of 1,000 samples/sec. Not all time points 
will be available, however, because of data transmission loss from the SRB-DTV. (Data that are not good, 
such as those that go beyond the calibration range or that result from a synchronization problem, will be 
rejected and data from the next time point will be read. Also, if M04, M05, or M06 have changed more 
than 5 deg/sec from the previous time point or if M56, M67, or M69 have changed more than 3 kip from the 
previous time point, then those data will be thrown out and the next time point of data will be read.) 

The data are then run through a fourth-order digital Butterworth low-pass filter whose desired break frequency 
is specified on the SRBNAM name list card. Subroutine COEF is used to compute coefficients for the filter. 
The values of PDOT, QDOT, and RDOT (p, q, and f) in deg/sec 2 -g are then computed between two con- 
secutive time points from filtered M04, M05, and M06 values; PDOT, QDOT, and RDOT are often filtered 
themselves. Accelerometer corrections are then applied to the filtered values of the three accelerometers 
(MOl, M02, M03) as follows: 

MOlc = M01/cos(2 .75 x 0.017452) 

M02c = M02 - 149 . 12 RDT 

M03c = [M03 - 149 .12PDOT - 2 .0 QDOT] /cos (2 .75 x 0 .017452) 

The seven hook loads Ri are then computed as shown in table C-l. Finally, the eight stresses a,- arc computed 
as shown in table C-2. (Filtered M56, M67, and M69 values, however, were substituted for Va, Vbl< and Vbr 
in the equation.) Time in total seconds and the eight stresses arc then stored on the permanent file. Another 
time point of data is then processed, and the process continues until the interval has been completed. 

The job card setup for running SRBHOOK is as follows: 

Jobname, Tseconds. 

USER (SRB1,FDPS) 

CHARGE (45,37) 

DEFINE (TAPE4=dfpf n/CT=S, M=R) 

dfpfn = pfn of data file of stresses to be created by SRBHOOK 
SWITCH (3) 

GETPF (LGO=SRBHOOK) 

GETPF (PROCFIL/UN-CSDO) 

BEGIN (TASK, , LIST, pname) 

pname = procedure name of the desired FLIDAB 


7/8/9 card (end-of-record card) 

SRBlnnn 

nnn = SRB flight number of desired FLIDAB 

$ SRBNAM $ 

$ TABLE 1 $ 

STABLE 2 $ 

6 / 7 / 8 / 9 card (end-of-input card) 
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SRBHOOK 


Revised for Flights 8 and 9 

As previously mentioned, during the entire SRB-DTV flight-test program there were only three channels avail- 
able for recording the strain gage outputs. From flight 1 up to flight 7 the three available channels were used to 
measure Va, Vbl> and Vbr (Sa, Dcl, Dcr, and Scl were not measured). In the original SRBHOOK program the 
measured hook loads Va, Vrl, and Vbr were used directly to compute some of the stresses. The nonmeasured hook 
loads Sa , Dcl, Dcr, and Scl > however, were approximated using accelerometer data, and then the other stresses 
were computed. 

For flights 8 and 9, the two channels used for measuring Vbl and Vbr were switched to measure Dcl and 
Scl instead. These measured values of Dcl and Scl were found to be significantly different from the values 
approximated from accelerometer data on flights 1 to 7. It was obvious that the accelerometer outputs were just 
too poor and unreliable to use for approximating hook loads. Therefore, it was decided to revise the SRBHOOK 
program for flights 8 and 9 to read only M56 (forward hook vertical load), DSSCL (left rear hook side load), and 
DSDCL (left rear hook drag load) and compute the stresses that were a function of them: 


<ti = 7.3522 M56 

<74 = 0.9688 ( —0 .0545 14 DSDCL x DSDCL + 6 .3481 DSDCL) 

+ 0.2479 (-0.3436 DSSCL x DSSCL -9 .9653 DSSCL) 
a 5 = 2.18675 DSDCL 

<76 =7.1420 (0.9688 DSDCL -0.2479 DSSCL) 

07 = 7.3055 DSDCL 
erg =4.2141 DSDCL 

A list of the revised SRBHOOK program is included in the following pages. 

THIN Program 

I. Functions 

A. Reduces the data sampling rate from 1 ,000 samples/sec to any specified number per second (50 sam- 
ples/sec is currently used). 

B. Interpolates new values where needed. 

C. Provides capability to write file header records for later identification, verification, or documentation. 

D. Provides a listing of all data dropouts, their duration, first and last values, and sums of the time loss in 
two categories: 

1. Total time loss 

2. Total of time losses greater than 1 sec. 

E. Echoes input control cards and file parameter identifications and comments. 

F. Lists any requested amount of data from the input file. 

G. Lists any requested amount of data as they go to the output file. 

II. Program organization 
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A. After declarations and initialization, the control cards are read and written out to the printer. 

B. The header cards are then read, written on TAPE3, and printed. The page headings are written to the 
print files. 

C. To separate data for the three listings (input from TAPE4, output to TAPE3, and dropouts), three discrete 
print files arc created: (1) the default printer for INPUT, (2) TAPE6 for output, and (3) TAPE7 for the 
dropout list. Appropriate headers arc written on the three files at this point. 

D. Because of the requirement to have two sets of data in memory whenever interpolation is required, the 
processing is done with two nearly identical sections of code working in a flip-flop mode and a special 
first data record section to initialize the first locations. 

E. In each of these sections TX (time expected) is compared with T. In this record T stands for the time 
found on the last data record read in and may be eitherTl orT2, depending on which section is processing 
the record. When TX is greater than T, another record is read and tested. If TX is equal to T, the record 
associated with T is written to the output tape and TX is increased by DTX (the desired delta time 
between output records). When TX is less than T, new values are interpolated for each parameter at 
each value of TX + N x DTX that occurs between T and T - 1 . The value TX is then incremented by 
DTX. During this procedure tests are being made to determine whether the input or output records, or 
both, are to be printed and when the limit time or end of file has been reached. 

F. On end of file or time limit, appropriate comments are written to the print files along with the final values 
and the summations. 


III. Input 


A. TAPE4 (created by the SRBHOOK program) 


B. Control cards (NAMELIST format) 


$CCDS 

NPARMS 

NCRDSP 

NCMTS 

NSKIP 

TSTART 

TSTOP 

DTX 

PTLIM 

PNLIM 

NDLIM 

NSIGMA 

$END 

Parameter 1 ID 
Parameter 2 ID 


Identifies the group 

Integer number of parameters on file 

Integer number of identification lines per parameter 

Integer number of comment lines 

Integer number of unused parameters 

Real time (sec.) to start processing 

Real time (sec.) to stop processing 

Real time (sec.) between output data parameters 

Real limit time for printing input data 

Integer limit count for printing input data 

Integer limit count for printing output data 

Integer number of parameters on file 

Group termination 

1 - NCRDSP identifying a parameter 
1 - NCRDSP identifying a parameter 


NPARMS 

Comments 1 


NCMTS comment cards 
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IV. Output 

A. FILE (TAPE3) 

1. Header format 


Record 

Column 

Mnemonic 

Type 

Description 

1 

1 to 80 

TITLE 

Character 

— 

2 

1 to 10 

NCHANS 

Integer 

Number of data parameters 

2 

11 to 20 

NCRDSP 

Integer 

Number of lines/parameter 

2 

21 to 30 

NCMTS 

Integer 

Number of comment lines 

3 

1 to 80 

CARD 

Character 

NCHANS sets of NCRDSP lines 





each describing the data 





channels 

2+NCRDS x NCHANS 



3+NCRDS x NCHANS 



2+NCRDS x NCHANS 



+NCMTS 





1 to 80 

CARD 

Character 

Comments 


2. Data format 


Record 

Heading 

1 

TIME 

2 

SIGMA 1 

3 

SIGMA 2 

9 

SIGMA 8 


B. LISTING 

1. Echo of the card input data 

2. Header records from TAPE9 

3. First PNLIM points 

4. Last record from TAPE4 

5. Last record output to TAPE3 

6. Postmortem dump of all variables 

7. First NDLIM points going to TAPE3 

8. TIME and sigmas spanning data dropouts 

9. Sum of all dropout periods 

10. Sum of dropout periods greater than 1 sec 
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V. Job setup 


THIN, T10000 . 

USER (user name, user ID) 
CHARGE (charge numbers) 
MAP, ON. 

FTN5 (ET, LO, DB=PMD) 
NOEXIT. 

ATTACH (TAPE4=input file 
name) 

LGO . 

REPPF (TAPE3= output file 
name) 

NOEXIT . 

REWIND , TAPE 6 , TAPE7 . 

COP YCF, TAPE 6, OUTPUT . 
COPYCF, TAPE 7 , OUTPUT . 
7/8/9 

PROGRAM DECK 
7/8/9 

CONTROL CARDS 
HEADER DATA CARDS 
6/7/8/9 


Job name and time estimate 


Optional 

Compile with cross 
references 

Prevents abort on forced 
dump 

Gets input file 

Starts program operation 
Names and catalogs the 
output 

Prevents an error condition 
from aborting the following 
commands 

Commands to move the 
temporary output files to 
the print file 
Separator 

Separator 

End-of-job indicator 


CONCAT Program 

I. Functions 

A. Concatenates files containing separate parts of a flight into a single file. 

B. Places all the header records from the various files onto the front of the first file on the output tape. 

C. Will remove any number of “bad” initial values that get through the filter algorithm during the startup 
process. 

D. Can be made to print any distance into each of the files for verification. 

II. Program organization 

A. The PARAMETER statement is used to assign values to the constants which control the program’s op- 
eration. 

B. CHARACTER strings and REAL arrays are defined in the next group of statements. 

C. The loop beginning with DO 35 reads the header records from each input file, prints them, and then 
writes them to the output file, TAPE10. TAPE9 is used as an intermediate transformation medium. 
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D. The section between DO 37 and 37 CONTINUE eliminates the first FP data values from the first file 
(that is, the first flight segment). 

E. DO 45 through 45 CONTINUE is the main data copy algorithm and logically includes the section from 
statement number 62 to the end of the program. 

1. In the DO 45 statement N is progressively incremented by 1 from 11 through 10 + NTS and specifies 
the file being read from. 

2. Statements 38 through 39 read from tape unit N, increment the input and output counters (NIN and 
NOUT), print (if the number of inputs is within the range specified by the print control parameter 
NPP), and write the data out to TAPE 10. 

3. On sensing an end of file, processing commences at statement 62 where the last data scan is printed, 
followed by the title line for the next file. 

4. The loop DO 64 through 64 CONTINUE reads and prints but does not save the required deletion 
records. 

5. Statement 65 through GO TO 65 provides a fill-in for any time gaps occurring between consecutive 
files by copying the last data scan over and over with an incremental time value until the new time 
is reached. 

6. Transfer then goes to the beginning to the DO 45 loop, which copies the rest of the file to TAPE10. 

F. Statement numbers 51 and 52 write appropriate messages when input errors or unexpected end-of-file 
characters are encountered. 


III. Input 

A. Tape files 

The file containing the first portion of the flight is assigned to TAPE11 , the second segment of the flight 
is assigned to TAPE 12, and so on, for as many as needed. 

B. Program controls 

Four integer parameters must be set in the program to provide the necessary processing information. 

L Number of data channels (8) 

NTS Number of files to be combined (2 to ?) 

NFP Number of the first usable data scan from each file 
(that is, eliminate NFP-1 scans) 

NPP Number of points to be printed from the beginning of each input file 

IV. Output 

A. A FILE denoted as TAPE10 is created that contains all the headers and data from the input files except 
those data scans specifically removed from the beginning of each input file. It also contains the “fill” 
data between the end of each original file and the beginning of the next. 

B. A LISTING is produced containing all the header records, the first NPP scans from each input file, the 
last scan from each file, and notification of read errors or premature end-of-file statements detected. 
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V. Job setup 


CONCAT, T10000 . 

USER (SRB1 , FDPS ) 

CHARGE (45, 37) 

FTN5 (ET, LO, DB=PMD) 

GETPF (TAPE11=S0705T1) 

GETPF (TAPE12=S0705T2) 

GETPF (TAPE13=S0705T3) 

GETPF (TAPE10+N=S0705TN) 

LGO . 

CATLIST. 

REPPF (TAPE10=S0705CX/CT=S) 
CATLIST. 

7/8/9 

PROGRAM DECK 

S/l/8/9 


Job name and time estimate 
User name and user 
identification 
Charge numbers 
Requests a FORTRAN compile 
Gets first part of flight 
Gets second part of flight 
Gets third part of flight 

Gets Nth part of flight 
Starts program execution 
Requests listing of file 
names 

Catalogs the newly created 
file 

Requests listing of file 
names (for comparison) 
Separator 

END-OF-JOB 


CCLAMP Program 

(Compute crack length at maximum points) 


I. Functions 

A. Computes crack length at every maximum and minimum value detected on each of the requested param- 
eters (DIRECT method). 

B. Extrapolates crack growth each time a string of invalid or missing data points is encountered. 

C. Uses data from other flights to compute crack growth (INDIRECT method). 

D. Computes remaining number of flights for each part. 

E. Compiles statistics on flight profile. 

F. Provides results of intermediate calculations and supporting data for validity checking. 

II. Method 

Data scans are read from the input file and checked against the various time boundaries. If the associated time 
is within the limits, processing continues. 

When a data field containing either a maximum or a minimum is detected on a parameter, it is held until the 
next extreme point in the opposite direction occurs on the same parameter. Then the maximum and minimum 
are passed to the CRACK subroutine. When the next extreme on the same parameter occurs, the process is 
repeated. 

When the maximum occurs earlier in time than the minimum, the slope indicator is set to — 1 ; in the reverse 
order it is +1. 
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On each subroutine call for a half-cycle calculation, the crack length multiplier is set to 0.5. 

If a time lapse greater than a specified limit (0.1 sec) occurs, extrapolation is performed by extending the 
average slope from the previous check point. 

Three types of time boundaries determine (1) the starting-stopping time, (2) segment limits for extrapolation 
and intermediate output, and (3) duration between plotting points. 

A. Peak detection 

A storage location, SIGN(J), is set aside for the sign of the slope of each parameter as it exists on the 
current record in memory. The value J indicates the parameter number. SIGN(J) is +1 when the slope 
of parameter J is positive and — 1 when the slope is negative. Another set of locations, SAVL(J) is set 
aside to hold the previous value of each parameter J. 

1. Detection of maximum occurs when the slope SIGN(J) is positive (+1) and the current value of 
the parameter FLIM(J) is less than the previous value SAVL(J). The value in SAVL(J) is used for 
calculation and is then replaced with the subsequent value in FLIM(J). The value of SIGN(J) is 
changed to —1 and the peak counter N(J) is incremented. 

2. Detection of minimum occurs when the slope SIGN(J) is negative (—1) and the current value of 
the parameter FLIM(J) is greater than the previous value. The previous value SAVL(J) is used for 
calculation then replaced with the current value FLIM(J). The variable SIGN(J) is then set to a 
+1 value. 

B. Subroutine CRACK is used to compute crack growth H and a new crack length CCL for each half-cycle 
by the following method. (Symbols in parentheses reflect mathematical variables used elsewhere in this 
report.) 

1 . H = CRx DELTK< NB AR ) x OMR*- SHX > x 0 .5 ( da/dN ) 

where: 

CR is the cycle ratio, 

NB AR is the slope of log (da/ dN ) as a function of log (DELTA K), and 

SHX is the strain hardening exponent. 

2. DELTK = DK1 x (SMAX - SMIN) x V / CCL 


where: 



SMAX 

is the current maximum. 

( &max) 

SMIN 

is the current minimum, 

( ^rm'n) 

CCL 

is the current crack length, 

(a) 

DK1 

is a constant computed once at the beginning of the program with values that do not 


change during the run, 


DK1 

_ CLP x FMF x \/pT 
\/CSF ~ 


CLP 

is the crack location parameter. 

(A) 

FMF 

is the flaw magnification factor, 

(M fc ) 

CSF 

is the crack shape factor, and 

(Q) 

PI 

is the constant 3.14159265359. 


OMR = 

1 - SMIN/SMAX 

(1 O’min/ CT mo x) 

Several 

nonvalid conditions are tested for and eliminated. 

These occur when OMR is negative 


or when DELTK or SMAX is zero or less. NC(J) is incremented by one each time any of these 


(m) 

(n) 

(A K) 


58 



conditions occur on parameter J. Printing of the calculation results and many of the input values 
occurs when the peak count N is less than the print control variable NCTL and whenever one of 
these nonvalid conditions exists. 

5. Some variables are passed through to the subroutine that do aid in interpreting the information being 
printed although they are not needed for calculations. These variables are: 

TMIN the time at which the minimum occurred, 

TMAX the time at which the maximum occurred, 

M the slope direction indicator, and 

NIM the number of the last record read prior to calculation. 

The value of L is computed to direct the print to a separate file for each parameter. 

C. Timing 

1 . START and STOP indicate the time when processing is to start and to stop. Zero values imply that 
the first time and the last time (respectively) in the file are to be used for these values. If the start 
time is earlier or the stop time is later, or both, than any file data, the crack length for each sigma is 
extrapolated appropriately. 

2. End-of-segment (EOS) times force the printing of intermediate results and provide reference points 
from which extrapolation occurs. 

3. Delta time (DELT) is an incrementing value for the TTIM, which forces the output of crack lengths 
to the plotting file. 

D. Extrapolation 

Although some of the short-duration data losses have little effect, others can exist for several seconds, 
and required compensation must be made for missing information. 

Extrapolation with respect to time is performed based on the average growth rate from the previous 
checkpoint. Checkpoints may be specified each time the type of flight activity changes, such as start 
taxi, stop taxi, resume taxi, begin takeoff run, airborne, or reduce to cruise power. The default start 
of interpolation occurs when more than five consecutive second-derivative values are less than epsilon 
(EPSL) (currently 0.1 x 10 _5 ) on the specified data channel NOTC (currently channel number 1) (see 
fig. 21). 

“Back” extrapolation is performed at the start of some flights when the first recorded data occur after the 
first motion of the aircraft is recorded in the flight log. Back extrapolation is a linear extrapolation back 
to the start from the first checkpoint through the first recorded data point (see fig. 22). 

III. Program organization 

A. Main program 

Initial coding was written to permit a variable number of parameters (as), but as development progressed, 
the number of as was frozen at eight. In the interests of time and efficiency, some code (primarily print 
statements) was written for the fixed number. 

1. After declarations and initialization, the control cards are read and written out to the printer. Ad- 
ditional initialization is performed in the block containing statements 1 and 2, which needs to be 
reset after back extrapolation. DK1, computed in this block, is the portion of the DELTK that uses 
variable values not changed during the flight. The temporary print files are initialized with the 
appropriate header information in the block terminated with statement 4. 

2. TAPE4 (the file generated by CONCAT) is read for its header information, which is printed for 
verification, then the dropout list file header is written. 

3. The first data record is read from TAPE4 by a special section of code, beginning with statement 
15, which initializes the START and TTIM variables with its time value and initializes the SAVL 
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locations with its parameter values. This section causes the first part of the file to be bypassed when 
start time is later than the first times in the file, and it initializes the variables required for back 
extrapolation when the start time is earlier than the first time in the data file. 

4. The block of code containing statements 22 and 24 is an end-of-file procedure to update the printer 
variable and to compute the final slopes if the indirect method is being used. 

5. At statement 26, a TAPE4 read is again performed with the time going into a variable named TIME 
and the parameter values going into an array named FLIM. Since TIME should be greater than 
TTIM, the crack sizes (initial values at this point) will be printed and written to files named TAPE6, 
TAPE3, TAPE7, and TAPE8 for listing and for later use by the CPLOT program. 

6. Just ahead of statement 28, TIME is compared with the current delta time terminator TTIM and 
the current end-of-segment marker (EOS). When either of these or an end-of-tape condition exists, 
indirect computed crack lengths CCLH are formed for those channels having a positive SLOPE 
value. Time and all the CCLH values (both direct and indirect) are written to plot (TAPE3) and 
print (TAPE6) files under control of format statements 101 and 100. 

7. The statement DO 40, J = 1 , NSIGMA starts a loop that analyzes each of the parameters (SIGMAS) 
one at a time. 

8. If the SLOPE value corresponding to the current SIGMA and flight segment is a positive value, 
direct calculation of crack size is bypassed, and the DO control integer J is incremented. When the 
SLOPE value is zero or negative (a nonpermissible physical condition), calculation proceeds down 
through statement 40. 

9. The variable SIGN is checked for slope; because it has had no chance to be changed from its ini- 
tialization to zero, the tests for greater than zero and less than zero will fail, and computation will 
continue at the test for SIGN = 0. Here any change in amplitude of each parameter is detected and 
is used to establish the initial sign of the slope and either the maximum or the minimum test value 
for the next record. 

10. After each parameter has been checked, control is transferred back to the data read statement (state- 
ment 26). When SIGN ^ 0, the comparisons between FLIM and SAVL become significant. When 
SIGN is positive and the last data point (FLIM) is less than the previous point (SAVL), a maxi- 
mum has been passed. The time and data in storage from the previous read are input to subroutine 
CRACK, which is called in the positive half-cycle mode. SIGN is set to a value of — 1 . 

11. When SIGN is negative and the change in value of a parameter is positive, a minimum has been 
passed. SIGN is set to a +1 , the data from the previous read are input to the subroutine CRACK, 
and it is called in the negative half-cycle mode. 

12. After an initial slope direction has been established, all points not qualifying for maximum or min- 
imum are ignored except for the test for invalid data and dropouts. 

13. Because the CONCAT program filled missing data areas with linear interpolated values, these ar- 
eas are detected by testing their second derivatives for values less than a very small EPSL. When 
the number of consecutive dropouts is greater than a specified number (NEPCS.GT.NEPS), linear 
extrapolation of crack size is computed by extending the slope from the last checkpoint through the 
last computed point. NEPS is currently set to 5, which is 0.1 Hz at the 5-cycle cutoff frequency and 
50 samples/sec. 

14. Transfer then goes to the top of the loop (statement DO 40) for analysis of the next parameter. 

15. After each parameter has been processed (as described in paragraphs 7 to 14) control goes back to 
the statement labeled 26 to bring in another record from TAPE4. The input record is printed if the 
record count is not greater than the limit LPRINT. The current crack size computed from each data 
parameter is output if the new time TIME is equal to or greater than the last value TTIME. The 


peak detection process is repeated for each channel as before. This process continues to the end of 
TAPE4 file. 

16. An end of file, an exceeded STOP time, or a zero end-of-segment (EOS) value cause all the outputs 
to take place, and control is moved to the statement labeled 50 where summations are formed and 
printed. Calls to subroutine HMS convert the times from seconds to hours, minutes, and seconds. 

17. One more linear interpolation is performed on the crack sizes if the last EOS time is greater than 
the last time in the data file. 

18. Based on the peak values from this flight, the remaining number of flights is calculated for each 
parameter from the crack lengths. 

19. Peaks from previous flights are compared, and remaining flights are calculated again based on the 
peak values detected on this and previous flights. 

20. The MOD array is used to modify the selection of slopes that are output for use in a subsequent flight 
when the indirect method is required but the flight pattern (that is, starts, stops, and taxi periods) 
does not match the current flight. 

2 1 . An output listing of all values to be entered in the subsequent flight is produced along with a relisting 
of the slope values with reduced precision but in more readable format. 

22. An input-selectable postmortem dump and FORMATS constitute the final statements in the main 
body of the C CLAMP program. 

B. Subroutine CRACK 

1. Initialization 

a. PARAMETER LOC size is set to 8. 

b. REAL arrays, including NBAR, which would normally be INTEGER mode, are declared to be 
size LOC specified in the PARAMETER statement 

c. INTEGER arrays are declared to size LOC, FLAG is declared LOGICAL and set to FALSE. 

d. The COMMON array’s order and size is declared consistent with the MAIN program’s COM- 
MON array. 

e. Output file control variable L is set to 10 plus the a number J for which crack length is being 
computed. 

2. Delta K (DELTAK) is computed using SMAX (the peak amplitude), CCL (the previously computed 
crack length for this location), and DK1 (a constant formed from input variables that do not change 
during the flight). 

3. OMR (1 — R) is computed using a minimum and a maximum (SMIN and SMAX). 

4. The various error conditions are checked, and if any exist, the growth increment H is set to zero, 
the FLAG is set TRUE, the error counter for this channel is incremented, and control is transferred 
around any further calculations to the print output. 

5. When the error conditions do not exist, H is computed and added to the total crack length. 

6. Input to and results of the calculation are written to an output file (TAPEL) for later printing when 
any of the following conditions exist: 

a. The error FLAG is TRUE. 

b. The number of printed lines requested (NCTL) has not been exceeded. 

c. A print control request (NPL) transmitted from MAIN is turned on (that is, equals 1). 

C. Subroutine HMS is a very short (three active lines) algorithm that converts time in seconds to hours, 

minutes, and seconds. 
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IV. Input 


A. Time and associated stress amplitude (SIGMA) data are entered from the tape format file created by 
CONCAT, assigned to local file name TAPE4. 

B. Constants and controls are entered from cards via the two NAMELIST groups TABLEX and TABLEC. 

1. TABLEX contains the values that are most frequently changed. This group is set up as follows: 


Name 

Type 0 

Content 

STABLEX 



ID 

C 

string of 1 to 70 characters in quotes 

NFLT 

I 

flight number 

NPRT 

I 

number of previous flight 

LOCNAM 

CA 

8 to 20 character strings, each in quotes 

LOCNR 

I 

sigma numbers for nonsequential assignment 

START 

R 

time to start process (default is start of file) 

STOP 

R 

time to stop process (default is end of file) 

DELT 

R 

time between crack size outputs 

NSKIP 

I 

spurious data channels to be skipped 

NSIGMA 

N 

number of as (data channels) 

LPRINT 

I 

number of input records to be printed 

NCTL 

I 

number of calculation lines to be printed 

IPMD 

I 

turns off the postmortem dump when IPMD = 0 

NOTC 

I 

number of test channel 

NEPS 

I 

dropout number threshold for interpolation 

EPSL 

R 

epsilon, minimum permitted value of second derivative 

EOS 

RA 

end-of- segment times (up to 10) 

CCLQ 

RA 

calculated crack length from previous flight 

SID 

CA 

segment identification 

PKVAL 

RA 

peak values (for each a) from prior flights 

PKTIME 

RA 

time in the flight when peak occurred 

NFLTPK 

IA 

number of the flight contributing the peak 

PKFLST 

RA 

peak flight start time 

NFLTP 

I 

number of flight previous 

NCYO 

I 

third code digit in plot and listing file 

IGNO 

I 

fourth code digit in plot and listing file 

NOFNCK 

I 

no flight number sequence check 

STARTP 

R 

start time for printing input data 

STOPP 

R 

stop time for printing input data 

SLOPE 

RA 

10 x 8 array of slopes for indirect calculations 

MOD 

SEND 

IA 

10 x 8 array to modify the slope’s sequence 


°C denotes CHARACTER, R denotes REAL, I denotes INTEGER, and A denotes ARRAY. 


2. TABLEC contains values that are less frequently modified. They are all REAL each having as many 
entries as there are parameters (currently 8) for each item name (except the group name $TABLEC 
and terminator SEND). 
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STABLEC 

CLP 

crack location parameter 

FMF 

flaw magnification factor 

CSF 

crack shape factor 

CR 

cycle ratio (from heat treating) 

NBAR 

stress intensity factor exponent 

SHX 

strain hardening exponent 

CAPC 

initial crack size from proof test 

SLIM 

sigma limit 

SEND 



V. Output 

A. Printer 

1. The control card input, TAB LEX and TABLEC, are echoed back for verification. 

2. The header records from TAPE4 are listed for their verification. 

3. Crack size summaries at the end of each flight segment are printed in order of occurrence. 

4. A summary of the entire flight lists flight duration, crack size growth rate, methods of calculation 
(direct calculation, interpolation, or indirect calculation), and number of flights remaining as deter- 
mined by various criteria. 

5. The next set of printed values is TABLEP, which is also the file punched out for entry into the next 
flight. 

6. TABLED shows the values of the remaining flights test criteria. 

7. The final directly printed table shows the crack growth rates for each cr for all flight segments. 

8. Several additional listings may be produced or eliminated under program input and control card 
setup by creating and then copying intermediate files to the printer. 

B. Intermediate files 

1. TAPE1 is actually the print file as described in the preceding section (V.A.). 

2. TAPE2 is used for temporary storage and for reformatting the header data on TAPE4. 

3. TAPE3 contains the crack sizes as functions of time, used primarily for plotting. 

4. TAPES is used to collect any error messages generated during calculations. 

5. TAPE6 contains crack sizes as functions of time, the same as TAPE3 except that TAPE6, being 
designed for the printer, contains a title block and column headings not included on TAPE3. 

6. TAPE9 contains a listing of the missing data areas (dropouts) for subsequent printing or plotting. 

7. TAPE11 to TAPE18 contain the input, output, and intermediate values associated with the crack 
calculation subroutine CRACK. Calculations of cti crack size are written to TAPE11, calculations 
of <72 are written to TAPE 12, and so on. 

8. TABLEP output in the form of punched cards is used to transmit final crack sizes, amplitude, and 
location of the highest peaks detected and crack growth rates (SLOPES) for use in the next flight. 
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Job name and time estimate 


VI. Job setup 

CCLAMP, T10000 . 

USER (username, Password) 

CHARGE (account number) 

MAP, ON. Provide a cross-reference 

map 

FTN5 (ET, LO, DB=PMD) Compile the source file 

* GETPF (TAPE4=f ilename) Get the file to be 

processed 

NOEXIT. Recover lists up to point 

of failure 

LIBLOAD (SRVLIB, GETJN) 

LGO. Start computing 

REWIND , TAP El , TAPE 5 , TAPE6, TAPE 7, 

TAPE8, TAPE9, TAPE10 . 

REWIND, TAPE11, TAPE12, TAPE13, TAPE14, TAPE15, TAPE16, TAPE17, TAPE18 . 
COPYEI, TAPE1, TAPE99 Copy information on the 

COPYEI, TAPE6, TAPE99 . first file named to the 

COPYEI, TAPE5, TAPE99 . second file named; the 

COPYEI, TAPE9, TAPE99 . range of data copied is 

COPYEI, TAPE11, TAPE99 . from the beginning to the 

COPYEI, TAPE12,TAPE99. end of information 

COPYEI, TAPE13, TAPE99 . 

COPYEI, TAPE14 , TAPE99 . 

COPYEI, TAPE15, TAPE99 . 

COPYEI, TAPE16, TAPE99 . 

COPYEI, TAPE17, TAPE99 . 

COPYEI, TAPE18, TAPE99 . 

PACK, TAPE99 . Removes file separators 

COPYBF,TAPE9 9, OUTPUT. OUTPUT is the print file 

REWIND, TAPE 9. 

* REPPF (TAPE9=S0705DL) Catalogs dropout lists for 

next program 

* REPPF (TAPE3=S0705CL) Catalogs crack lengths for 

next program 

7/8/9 Separator 

Program deck 

7/8/9 Separator 

Control cards 

6/7/8/9 End-of-job indicator 

^Indicates cards to be changed from one run to the next. 
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CPLOT Program 


I. Functions 

A. Generates a tape for use by CALCOMP plotters, model 925 plotter controller, to create annotated graphs 
(figs. 23 to 30). 

B. Draws any number of function-time curves on each sheet 

C. Creates multiple vertical and horizontal scales according to specified values or generates them automat- 
ically with calculations from data maxima and minima. 

D. Can draw a series of horizontal lines with vertical end markers to indicate such things as areas of inter- 
polated data. 

E. Can draw annotations at any angle or location. 

F. Can place a block of legend lines at any desired location. 

G. Will write a title at the top of each page. 

H. Accepts input data and control parameters from either cards or tape files in any combination. 

I. Prints five levels of diagnostic messages under control of an input parameter. 

II. Program organization 

A. PARAMETER, TYPE SPECIFICATION, and DATA statements set some nondefault variable types, 
array sizes, and default values. 

B. Control parameter input name lists 

1. CPDATA is the primary control variable input list. 

2. ZDATA is a secondary input list used only when more than one each horizontal or vertical scale are 
required. 

C. Call to TIME subroutine entry point TIME1 establishes a reference for measuring elapsed time periodi- 
cally throughout the processing so that excessively long procedures may be optimized in future runs and 
may provide a clue to input errors. 

D. The call to PLOTS and PLOT are required by the plot subroutines to initialize and establish the origin 
location; the call to FACTOR sets the unit distance increment. A default size of 2 cm is used in place 
of the 1-in. unit assumed in the CALCOMP writeup. All subsequent uses of the word “inches” actually 
references a 2-cm distance and are multiplied by 2/2.54 to obtain the true size; other spacing may be 
specified in the first set of control cards. 

E. Statements 1 to 10 store all control parameter sets on a file named TAPE1 for future recall. Messages 
and various levels of print echo are controlled by the IECHO parameter. TAPES is used for intermediate 
storage to permit the input to be printed exactly as entered, whereas subsequent input data displays show 
all the possible input parameters with their default values or carryovers from previous settings. 

F. The next block of statements down through number 25 load the data onto TAPE2 so that one independent 
variable may be paired with several dependent variables (that is; TIME as a function of SIGMA 1, TIME 
as a function of SIGMA2, . . . , TIME as a function of SIGMA8). 

G. Statements 29 to 32 control the positioning of TAPE2 looking for a data set whose FID (identification 
field) matches the one specified in the data set identification parameter XYCG. When no field identifi- 
cation is required (that is, when only one data set is in the file) XYCG may be set to “* 
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H. The next group of statements initialize the storage locations for the minimum and maximum X and Y 
values. Then beginning with statement DO 37 through 37 CONTINUE, each pair of XY coordinates is 
read and stored in the X and Y arrays. Input values may be adjusted by subtracting a constant (XMIN or 
YMIN) and then multiplying by another constant (XMPF and YMPF). When the MIN values are zero 
and the MPF values are 1.0, the input values are unchanged as stored in X and Y arrays. If the MIN 
parameters are set to the code value —9999 ., the first data point is subtracted from all others in the set. 
This capability permits the conversion of time in seconds from midnight to be plotted in elapsed time in 
minutes and crack length in inches to be shown as crack growth in microinches. The highest and lowest 
values of each X and Y array are saved in XHI, XLO, YHI, and YLO locations. Input is terminated in 
one of three ways: 

1 . an end of file detected on READ, 

2. the value of NPTS is exceeded, or 

3. die size of the X and Y arrays are exceeded (currently set at 1000 each). 

I. X scale factors specified by XSF are applied unless XSF is zero, in which case a scale factor is calculated 
from the highest (HI) and lowest (LO) X values saved during input. CALCOMP's subroutine SCALE 
does the actual calculation after the values are provided. 

J. Y scale factors are calculated in the same manner using YSF for the controlling parameter. 

K. The Y -axis length here becomes a test value since it cannot be zero if it is to exist at all. A YLEN = 0 
indicates that no axis, annotation, or titles are to be drawn on this pass. 

L. A TITLE of up to 80 characters is drawn at a fixed location — 10.5 “inches” up and starting 1 “inch” to 
the right of the Y axis. The letter height is 0.2667 “inches” (the inches in quotes are actually 2 cm). 

M. A vertical (Y axis) is drawn by line 43. The default location is at the left edge of the plotting area 
starting at location (0.0, 0.0). YXPG and YYPG control its position and may be set to other values. The 
characters stored in YNAME will be written vertically at the left of the scale. NYCH may be used to 
specify the number of characters (up to 40) in the string when positioning is to be controlled, otherwise 
the character string length (calculated by statements DO 42 to 42 CONTINUE) will be calculated to 
center the axis annotation. YLEN specifies the length of the Y axis in “pseudo” inches (2 -cm units) and 
may not be zero unless all scales and annotation are to be omitted. YIVAL is the annotation for the first 
tick mark of the axis. YSF(1) is the value increment between tick marks. 

N. When the parameter NZAX is greater than zero, NZAX auxiliary axes will be drawn. Z axes can be 
placed anywhere at any angle but the default orientation and location is vertical, 2 cm left of the previous 
axis. These may be changed by setting ZANGE, ZXPG, and ZYPG. The length and the scale factor 
defaults, ZLEN and ZSF, are the values used for the Y axis. ZIVAL, the first tick mark annotation, 
default is zero. ZNAME is used to enter an axis title up to 40 characters and NZCH is used as YNCH as 
in the previous section. 

O. The horizontal (X) axis is drawn starting at the default location of (0.0,0.0). The X-axis control param- 
eters XNAME, NXCH, XXPG, XYPG, XLEN, XIVAL, XSF, XMIN, XMPF, and LOCX have the same 
function as their Y -axis counterparts. 

P. Statements DO 54 through 54 CONTINUE control the generation of a LEGEND block of 1 to 20 lines, 
each up to 40 characters in length. 

Q. Statements DO 56 through 56 CONTINUE control the generation of 20 markers that may be placed at 
any location and at any angle. The placement (X and Y location), size, and angle of the 40 MARK 
characters is by the XYSAM variable array locations 1 to 4, respectively. 

R. Up to 20 vertical markers (VMARKs) are similar to the MARKs except that no angle is specified and the 
X location is specified in terms of the X value of the input data rather than page coordinates. XYSVM 
array locations 1 to 3 specify X value, Y location, and letter height, respectively. 


S. Up to 20 horizontal markers (HMARKs) are similar to VMARKs except that the Y location is specified 
in terms of the values on the Y axis. XYSHM array locations 1 to 3 specify X location, Y value, and 
letter height, respectively. 

T. Call to subroutine LINE produces a line plot of the pairs of data values in the X and Y arrays. The data 
points may be represented by centered symbols or connecting lines between points, or both. LTYPE 
and INTEQ control the line type and symbol selection. See CALCOMP line subroutine writeup for 
description of their use. Prior to calling LINE a test is made for “out of range” condition, which aborts 
the plot and prints a diagnostic message when the data fall more than two “inches” left or below (0.0, 0.0), 
more than 12 “inches” high, or more than 4 “inches” past the right end of the X axis. 

U. DROPOUT areas, where data were extrapolated, can be indicated on the plot. The statement 
IF(DOL.NE.0)THEN commences the section of code by testing for a nonzero value, which then in- 
dicates the height in pseudoinches above the X axis. 

V. Plot identification and tear marks for each plot are generated by the statements in the block between 
lines 80 and 8 1 and written vertically 4 “inches” to the right of the right end of the X axis. Plot se- 
quence, number, time, date, and JSN are computed or automatically picked up from the system. User 
identification may be entered as a four-character string assigned to the USER parameter. 

W. The reading of a new set of primary control parameters occurs in the program after statement 8 1 through 
the GO TO 29 statement and includes the following processes: 

1 . Location of the new page XNEW is increased to a point 8 “inches” past the end of the X axis 
(XLEN), and page and plot counters are incremented. 

2. Name list CPDATA is read and checked for end of file and read ERRORS. 

3. The ID file FIDA is searched for, and pointer K1 is set, on a match with the current field ID(FID). 

4. The ID file is again searched for, and pointer K2 is set, on a match with the data set identification 
parameter XYCG. 

5. Comparison between K1 and K2 determines the need of a rewind to locate the data requested in the 
new CPDATA parameter set. 

6. If the test of YLEN reveals it not to be zero, indicating an overplot, a call to PLOT reinitializes a 
new page. 

7. Processing then returns to statement 29 (see paragraph II.G) where the new data set is searched for 
and loaded into the X and Y arrays. 

X. Termination procedure 

1. Statement 90 reports on end-of-control-parameter-sets condition. 

2. Statement 95 provides notification of a tape error detected on the data file (TAPE2). 

3. Statement 96 moves the pen off of the page, flushing the buffer and writing a termination file (999) 
on the plot tape. 


III. Input 

A. Data files 

1. PRIMARY file format for TAPE3 is a file containing no header records having the following for- 
matted fields repeated in each record. 
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Program 

variable 

Field size 


name 

and type 

Application 

FID 

A5 

field for identification 

NCD 

15 

number of card (sequence number) 

NFPC(l) 

F10.3 

independent variable, X (time) 

NFPC(2) 

E15.9 

first dependent variable Y1 (cri ) 

NFPD(9) 

E15.9 

Eighth dependent variable Y8 (as) 


2. An optional auxiliary file, TAPE9, which may be attached and processed, contains the list of dropout 
areas. The format is X, 3I8.5F14.3, of which the first three integer fields and the last three floating- 
point fields are unused. The first two floating point fields contain the start time and the end time 
for a period during which extrapolation occurred. The file contains one record for each section of 
extrapolation during the flight. 

B. Control parameters (NAMELIST format) 

One or more sets of control parameters are required for each plot. The primary set is preceded by a 
SCPDATA field and followed by a SEND field. This set of controls can provide all the specifications 
for a single set of data points, the X and Y axes, the annotadon, and a dropout line. Secondary sets are 
used to specify additional axes ( Z axes). Secondary sets are preceded by a SZDATA field and followed 
by the standard SEND field. When more than one set of data points is to be plotted on the same sheet 
previously generated, a new primary set of controls is used. All axes and annotation are suppressed by 
setting the Y -axis length parameter YLEN to 0.0. Most control parameters have default values and all 
are saved from one plot to the next so that only new or changed values need be entered. 


Mnemonic 

Type 0 

Use 

Default 

General parameters 

SCPDATA 

— 

specifies beginning of a primary data set 

— 

FACTR 

R 

adjusts size of plot 

0.7874016 

NSETS 

I 

number of data sets 

1 

NFPC 

I 

number of fields per card (input line) 

9 

IECHO 

I 

diagnostic print control 

1 

LTYPE 

I 

type of line connecting points 4 

0 

ENTEQ 

I 

plotting symbol 4 

0 

XYCG 

5CH 

data set identification 

Blank 

NZAX 

I 

number of additional ( Z ) axes 

0 

DOL 

R 

dropout line, “inches” from X-axis 

0 

DATA 

4CH 

input file source TAPE or CARD 

CARD 

USER 

4CH 

computer user identification 

Blank 

ANOT 

40CH 

annotation to go at end of plot 

Blank 

NPTS 

I 

number of points to be plotted (0 = all) 

0 

XL 

R 

left margin spacing, “inches” 

4.0 

XT 

R 

not implemented 

4.0 

SEND 

— 

specifies end of each data set 

— 



Mnemonic 

Type 0 

Use 

Default 

Annotation parameters 

TITLE 

80CH 

graph title 

Blank 

LEGEND(l) 

40CH 

legend line 1 

Blank 

LEGEND(20) 

40CH 

legend line 20 

Blank 

LEGSYM(l) 

I 

legend 1, first char symbol indicator^ 

0 

LEGSYM(20) 

I 

legend 20, first char symbol indicator 

0 

NLEGS 

I 

number of legend lines 

0 

XYSLEG(l) 

R 

X location of first legend 

0.0 

XYSLEG(2) 

R 

Y location of first legend 

0.0 

XYSLEG(3) 

R 

size of legend letters 

0.0 

NHM 

I 

number of horizontal markers 

0 

HMARK(l) 

40CH 

horizontal annotation marker 1 

Blank 

HMARK(20) 

40CH 

horizontal annotation marker 20 

Blank 

XYSHM(1,1) 

R 

X location of first horizontal marker 

0.0 

XYSHM(2,1) 

R 

Y value of first horizontal marker 

0.0 

XYSHM(3,1) 

R 

size (letter height) of first horizontal marker 

0.0 

XYSHM(1,2) 

R 

X location of second horizontal marker 

0.0 

XYSHM(3,20) 

R 

size of 20th horizontal marker 

0.0 

NVM 

I 

number of vertical markers 

0 

VMARK(l) 

40CH 

vertical annotation marker 1 

Blank 

VMARK(20) 

40CH 

vertical annotation marker 20 

Blank 

XYSVM(1,1) 

R 

X value of first vertical marker 

0.0 

XYSVM(2,1) 

R 

Y location of first vertical marker 

0.0 

XYSVM(3,1) 

R 

size of first vertical marker 

0.0 

XYSVM(1,2) 

R 

X value of second vertical marker 

0.0 

XYSVM(3,20) 

R 

size of 20th vertical marker 

0.0 

NAM 

I 

number of annotation markers 

0 

MARK(l) 

40CH 

annotation marker, any angle 

Blank 

MARK(20) 

40CH 

annotation marker, any angle 

Blank 

XYSAM(1,1) 

R 

X location of first any-angle marker 

0.0 

XYSAM(2,1) 

R 

Y location of first any-angle marker 

0.0 

XYSAM(3,1) 

R 

size of first any-angle marker 

0.0 

XYSAM(4,1) 

R 

angle of first any-angle marker 

0.0 

XYSAM(1,2) 

R 

X location of second any-angle marker 

0.0 

XYSAM(4,20) 

R 

angle of 20th any-angle marker 

0.0 


Mnemonic 

TVpe° 

Use 

Default 

Horizontal axis parameters 

XNAME 

40CH 

X-axis title 

Blank 

NXCH 

I 

number of X axis title characters 

N 

XXPG 

R 

X-axis X starting location on page 

0.0 

XYPG 

R 

X-axis Y starting location on page 

0.0 

XLEN 

R 

X-axis length 

0.0 

XIVAL 

R 

X-axis initial value 

0.0 

XSF 

R 

X-axis value between tick marks 

0.0 

XMIN 

R 

X -axis minimum value (offset) 

0.0 

XMPF 

R 

X data multiplication factor 

1.0 

LOCX 

I 

location of field with X data 

1 

Vertical axis parameters 

YNAME 

40CH 

Y-axis title 

Blank 

NYCH 

I 

number of Y-axis title characters 

N 

YXPG 

R 

Y-axis X starting location on page 

0.00 

YYPG 

R 

Y-axis Y starting location on page 

0.00 

YLEN 

R 

Y-axis length 

10.00 

YIVAL 

R 

Y-axis initial value 

0.00 

YSF 

R 

Y-axis value between tick marks 

0.00 

YMIN 

R 

Y -axis minimum value (offset) 

0.00 

YMPF 

R 

Y data multiplication factor 

1.00 

LOCY 

i 

location of field with Y data 

2 

Auxiliary ( Z ) axis parameters 

SZDATA 

— 

specifies start of secondary data set 

— 

ZNAME 

40CH 

Z-axis title 

Blank 

NZCH 

I 

number of Z-axis title characters 

N 

ZXPG 

R 

Z-axis X starting location on page 

1.00 

ZYPG 

R 

Z-axis Y starting location on page 

0.00 

ZANGL 

R 

Z-axis angle 

90.0 

ZLEN 

R 

Z-axis length 

YLEN 

ZIVAL 

R 

Z-axis initial value 

0.00 

ZSF 

R 

Z-axis value between tick marks 

YSF 


“Type codes 

no values are associated with these control words 

I integer data, numeric values without decimal points 

R real data, numeric values, must have a decimal point 

CH CHARACTER data, a string of characters bracketed by single quote 

marks. The preceding numbers indicate the maximum size of the charac- 
ter string. 

‘Default codes 

N value is determined from actual length of string 

Blank initial contents are all blank characters 
YLEN initial value is taken from previous YLEN setting 

Y SF initial value is taken from the previous Y SF setting 

“Refer to CalComp user’s manual (California Computer Products, Inc., 1976) 
for description of line types and plotting symbols. 


IV. Job setup 


CPLOT, T10000 . J 
USER (SRB1 , FDPS ) U 
CHARGE (45, 37) C 
FTN5 (ET, LO, PL) R 
GETPF (TAPE3=S0705CL/UN=SRB1) G 
GETPF (TAPE9=S0705DL/UN=SRB1 ) G 
GETPF (PENLIB/UN=LIBRARY) G 
LIBLOAD (SRVLIB, GETJN) G 
LDSET (LIB=PENLIB, PRESET=ZERO) I 
LOAD (LGO) L 
EXECUTE S 
GETPF (PROCFIL/UN=CSDO) G 
BEGIN (PROCFIL/UN=CSDO) S 
BEGIN, WRITAPE, , PL, , , TAPE4 , TIME+, 
ENDIF (DOTHIS ) 


Job name and time estimate 
User name and user ID 


Charge numbers 
Requests FORTRAN compile 
Gets primary data file 

Gets secondary data file (dropout list) 

Gets the plot subroutines 

Gets a system subroutine 

Initializes memory locations to zero 

Loads the FORTRAN object file 

Starts the program execution 

Gets plot procedure file 

Starts plot output procedure 


7/8/9 

PROGRAM DECK 
7/8/9 

CONTROL CARDS 
6 / 7 / 8 / 9 


Separator 


Separator 


Job termination card 


★Indicates cards to be changed from one run to the next. 
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Table C-l. Coefficients of h(K)k load equation for hook loads Ri = R ai i] x + R b> rj y + Re, T] z + Rd,P+ R ei 4 + Rj,r 
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Shaft — left sway brace 0 0 0 4.7559 x 10 ~ 

Shaft — right sway brace 0 0 2.18675 x 10 ~ 3 0 

Left drag pin 0 0 0 7.1420 x 10“ 

Right drag pin 0 0 7.3055 x 10 _3 0 

Right drag shaft support lower 0 0 4.2141 x 10 ~ 3 0 
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Figure C-l. FREDS overview. 
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Figure C-2. Crack length calculation. 
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